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types  of  impulsive  sounds.  The  first  series  compared  N-waves  over  a range  of 
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response  tradeoff  between  rise  time  and  overpressure,  and  duration  and  over- 
pressure was  measured.  For  equal  judged  loudness,  10  ms  rise  time  required 
8 dB  higher  overpressure  than  for  1 ms.  Duration  had  little  effect  in  the 
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range  100  to  200  ms  but  at-  250  ms  noticeably  enhanced  the  loudness.  These 
results  confirm  those  measured  by  Shepherd  and  Sutherland  (except  for  the  250  ms 
duration),  and  extend  the  measurements  down  to  0.22  ms  rise  time.  --'vyi  ■ 

The  second  series  tested  certain  "flat  top"  sonic  boom  signatures,  which 
according  to  current  theory  could  be  generated  by  special  very  long  SST 
aircraft  designed  for  minimized  sonic-boom;  these  were  compared  for  loudness 
with  a reference  N-wave  (pfj  = 0.5  psf,  1 ms  rise  time,  150  ms  duration).  The 
subjective  loudness  was  found  to  be  dominated  by  the  front  (or  rear)  shock 
ApSH  while  the  maximum  peak  amplitude  Ap^^^  of  the  inclined  "flat  top"  had  only 
a small  effect.  The  results  for  equal  loudness  were  well  fitted  by  an  empirical 
law  ApsH  + 0.11  ApMAx  Apfj.  This  shows  that  for  equal  loudness  the  peak 
amplitude  Apjii^x  of  flat  top  signature  is  substantially  higher  than  that  of 
the  N-wave;  thus  for  equal  amplitude  the  flat  top  signature  is  the  quieter. 

The  third  series  compared  filtered  N-wave  signatures,  using  a high-pass 
digital  filter  with  an  unfiltered  N-wave  signature  (1  psf,  1 ms  rise  time, 

150  ms  duration).  Two  cut-off  frequencies  were  used:  25  and  50  Hz.  The 
amplitude  differences  for  equal  loudness  were  very  slight:  less  than  0.6  dB 
at  most.  Thus  the  'infrasonic'  low  frequency  content  of  sonic  boom  N-waves  - 
although  it  dominates  the  spectrum  - has  no  significant  influence  on  the 
subjective  loudness.  Similar  tests  with  annoyance  as  a judgement  criterion  has 
showed  a tendency  to  increase  annoyance  for  filtered  N-waves,  but  again  this 
was  very  slight. 

In  the  last  test  series  the  tradeoff  between  overpressure  and  duration  was 
found  for  idealized  quarry  blast  signatures  composed  of  sequences  of  25  ms 
long  pulses  with  0.22  ms  rise  time.  The  range  of  durations  extended  from  25  to 
400  ms.  At  the  short  durations  the  loudness  increased  2 dB  for  each  doubling 
of  duration;  above  100  ms  the  increase  was  progressively  lower,  approaching  as 
an  asymptote  the  level  of  continuous  sound. 

The  results  in  each  series  were  compared  with  theoretical  predictions  by 
the  method  of  Johnson  and  Robinson.  All  but  the  long-duration  quarry  blast 
judgements  were  found  to  be  in  very  good  agreement  in  terms  of  relative  loud- 
ness levels.  With  an  ad  hoc  - but  physically  plausible  - modification  (includin' 
adjustment  of  the  critical  integration  time  of  the  ear)  the  predictive  method 
was  extended  to  encompass  the  long  duration  signals  as  well.  Thus  the  applic- 
ability of  the  method  has  been  demonstrated  for  other  types  of  transient 
sounds  than  the  N-wave;  and  the  extension  to  the  method  appears  to  bridge  the 
range  between  impulsive  and  continuous  sounds  of  similar  spectral  content. 
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Abstract 


OA  loudspeaker- driven  simulation  booth  with  extended  rise  time  capability 
(down  to  0.22  ms)  has  been  used  for  subjective  loudness  tests  of  sonic-boom  and 
other  types  of  impulsive  sounds.  The  first  series  coii5)axed  N- waves  over  a range 
of  0.22  to  10  ms  rise  time,  100  to  250  ms  duration  and  from  0.5  to  4 psf  (the 
latter  for  the  longer  rise  times)  (24  to  192  N/j*^)  peak  overpressure.  The 
response  tradeoff  between  rise  time  and  overpressure,  and  duration  and  over- 
pressure was  measured  | For  equal  judged  loudness,  10  ms  rise  time  required  8 dB 
higher  overpressure  tl^  for  1 ms.  Duration  had  little  effect  in  the  range  100 
to  200  ms  but  at  259^s  noticeably  enhanced  the  loudness.  These  results  confirm 
those  measured  by^- Shepherd  and  Sutherland  (except  for  the  250  ms  duration),  and 
extend  the  measurements  down  to  0.22  ms  rise  time. 

^he  second  series  tested  certain  "flat  top"  sonic  boom  signatinres, 
which  according  to  current  theory  could  be  generated  by  special  very  long  SST 
aircraft  designed  for  minimized  sonic-boom-,  these  were  ccanpared  for  loudness 
with  a reference  N-wave  (pjj  = 0.5  psf,  1 ms  rise  time,  150  ms  duration) . '■T?He 
subjective  loudness  was  fovind  to  be  dominated  by  the  front  (or  rear)  shock  Zpsjj  j 
while  the  maxim\im  peak  an5)litude  of  the  inclined  "flat  top"  had  only  a 

small  effect.  The  results  for  equal  loudness  were  well  fitted  by  an  empirical 
law  + 0.11  APmax  ~ ^N-  This  shows  that  for  equal  loudness  the  peak 

tude  of  the  flat  top  signature  is  sxibstantially  higher  than  that  of  the 

N-wave;  ^us  for  equal  amplitude  the  flat  top  signature  is  the  quieter. 

Vlhe  third  series  compared  filtered  N-wave  signatures,  using  a high- 
pass  digital  filter  with  an  unfiltered  N-wave  signature  (l  psf,  1 ms  rise  time, 

150  ms  duration).  “^Two  cut-off  frequencies  were  used:  25  and  50  Hz.  The  anpli- 
tude  differences  for  equal  loudness  were  very  slight:  less  than  0.6  dB  at  most. 
Thus  the  'infrasojiic'  low  frequency  content  of  sonic  boom  N- waves  - although  it 
dominates  the  spectrum  - has  no  significant  influence  on  the  subjective  loudness. 
Similar  tests  with  annoyance  as  a judgement  criterion  has  showed  a tendency  to 
increase  annoyance  for  filtered  N- waves,  but  again  this  was  very  slight. 

<^'In  the  last  test  series  the  tradeoff  between  overpressiare  and  dviration 
was  fotmd  for  idealized  qixarry  blast  signatures  composed  of  sequences  of  25  ms 
long  pulses  with  0.22  ms  rise  time.  The  range  of  durations  extended  from  25  to 
400  ms.  \At  the  short  durations  the  loudness  increased  2 dB  for  each  doubling 
of  diu-atmon;  above  100  ms  the  increase  was  progressively  lower,  approaching  as 
an  asyupwte  the  level  of  continuous  sound. 

^The  results  in  each  series  were  compared  with  theoretical  predictions 
by  the  method  of  Johnson  and  Robinson.  All  but  the  long-duration  quarry  blast 
judgements  were  foiaid  to  be  in  very  good  agreement  in  terms  of  relative  loudness 
levels.pl  With  an  ad  hoc  - but  physically  plausible  - modification  (including 
adjustment  of  the  critical  integration  time  of  the  ear)  the  predictive  method  was 
extended Y©  encompass  the  long  duration  signals  as  well.  Thus  the  applicability 
of  the  method  has  been  demonstrated  for  other  types  of  transient  sounds  than  the 
N-wave;  aiid  the  extension  to  the  method  appears  to  bridge  the  range  between 
impulsive  jmd  continuous  sounds  of  similar  spectral  content. 
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1.  INTRODUCTION 
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It  is  increasingly  evident  that  persistent  or  intense  noise  has  a 
complex  and  often  hannfvil  effect  on  human  beings.  This  is  especially  true  in 
the  case  of  sharp  transient  or  in^julsive  sounds.  Exaii5)les  are  the  noise  of 
Jack-hammers,  industrial  percussive  processes,  quarry  blast  waves,  g\mshots, 
and  sonic  booms.  Relevant  research  on  such  noises  has  been  both  ad  hoc  and 
scientific;  some  major  objectives  have  been  (i)  to  find  valid  methods  for 
predicting  the  loudness  of  a given  transient,  given  its  waveform  or  spectrum; 
(ii)  to  establish  damage-risk  criteria;  (iil)  where  practicable,  to  reduce 
the  anplitude  at  the  source,  or  to  modify  its  waveform  for  reduced  loudness 
(e.g.,  as  in  recent  sonic  boom  research). 

There  do  exist  several  proposed  methods  for  predicting  the  subjective 
loudness  of  sonic  boom  signatures  (e.g.,  [l],  [2]).  Being  based  on  specific 
spectra,  these  have  the  potential  for  application  to  arbitrary  in5)ulsive  wave 
forms.  The  aim  of  the  present  investigation  has  been  to  test,  and  revise  if 
necessary,  these  predictive  methods  in  application  to  a variety  of  inpulsive 
waveforms . 

A loudspeaker- driven  simulation  booth  ([3]- [5])  with  extended  rise- 
time capability  (down  to  0.22  ms)  has  been  lysed  for  subjective  loudness 
Judgements  for  can5)axison  with  the  theory.  A new  conputer- aided  scheme  was 
developed  to  predistort  the  electrical  input  signal  to  counteract  the  loud- 
speaker-booth  distortion  . The  essence  of  the  idea  is  to  alter  the  electrical 
input  spectriutt  by  the  inverse  of  the  loudspeaker-booth  transfer  function 
(cooplex  frequency  response),  thus  effecting  cancellation. 

Foirr  series  of  tests  were  made,  each  of  interest  in  its  own  right, 
aside  from  providing  a predictive  test.  Three  of  these  had  to  do  with  sonic 
boom  waveforms  - standard  and  nonstandard  - and  a foiu-th  with  simvilated  quarry 
blast  waves. 

The  first  series  deals  with  the  standard  sonic  boom,  with  pressiire- 
time  signature  like  a letter  "N".  The  steepness  of  the  front  and  reeir  of  the 
"N-wave"  (which  are  shocks^  are  described  by  their  'rise-time',  and  the 
interval  between  - length  of  the  N - by  the  'duration' . The  tests  attenpt  to 
check  earlier  enpirical  dependence  [6]  of  the  subjective  loudness  on  rise  time, 
duration,  and  N-wave  aaplitude. 


The  second  series  deals  with  certain  'flat-top'  sonic  boom  signatures: 
they  look  like  a letter  "N"  with  the  top  and  bottom  cut  off:  Accor- 

ding to  current  theory  ([7]-[l0])  such  'flat  top'  signatxares  wouldbe  generated 
by  a special  family  of  very  long  (^  300  ft)  SST  aircraft  designed  for  minimized 
sonic  boom. 

The  third  series  explores  the  contribution  of  the  'infrasonic'  low 
frequency  content  of  sonic  boom  N- waves  to  their  subjective  loudness,  or, 
alternatively,  annoyance.  For  the  test  N-waves  the  frequencies  below  25  Hz 
and  below  50  Hz,  respectively,  axe  cut  off  by  digital  filters  simulating  sinple 
RC  circuits.  Conparison  is  made  with  a standard  unfiltered  N-wave. 
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The  fourth  series  deals  with  the  loudness  of  simulated  quarry  blast 
waves.  These  axe  generated  as  multiple  "shock"  seqviences  (simulating  echoes), 
with  e^qponential  decay  between  shocks:  The  length 
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(number  of  shock-decay  repetitions)  of  the  sequences  is  varied  over  a wide 
range:  from  short  iiqpulsive  bursts  to  almost  continuous  sound. 

For  each  of  these  test  series,  using  human  subjects  for  loudness 
Judgements,  the  results  were  con^jared  with  loudness  predictions.  For  the 
first  series  both  the  Zepler-Harel  [1]  and  Johnson-Robinson  [2]  predictive 
methods  were  used.  The  latter  appeared  to  have  somewhat  stq)erior  accuracy 
over  the  full  range  of  rise  times,  as  well  as  greater  sin^jlicity.  Thus  only 
the  Johnson-Robinson  method  was  applied  in  all  the  later  work.  For  the  fourth 
series  - the  simulated  blast  waves  - an  ad  hoc  modification  was  made  to  the 
method  to  bridge  the  transition  from  in^juisive  to  continuous  sound. 
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2. 


REVIEW  AND  BACKGROUND 


2 .1  Introduction  - Effects  of  Noise  on  People 

A commonly  used  definition  states  that  noise  is  mwanted  or  lanpleasant 
sound.  Thus  noise  is  another  pollutant  in  the  human  environment.  It  can  cause 
auditory  problems  by  masking  speech;  fiirther,  it  can  produce  "extra-auditory" 
effects  by  disturbing  physiologic  functions.  It  can  also  change  task  performance, 
disriQit  rest,  relaxation,  and  sleep.  People  exposed  for  long  periods  of  time  are 
can5)laining  more  often  that  noise-induced  stress  has  a degrading  effect  on  their 
health.  The  reported  problems  (e.g.,  Welch  and  Welch  [12])  include  neurologic, 
digestive,  and  metabolic  disorders,  cardiovascular  problems,  and  even  mental 
problems,  hypertension,  nervousness,  etc. 

It  is  well  known  that  high  intensity  sound  may  produce  temporary  or 
even  permanent  losses  of  hearing.  This  is  docimented  by  a large  number  of 
cited  examples  and  ptiblished  investigations.  The  main  task  of  those  investiga- 
tions was  the  production  of  a damage  risk  criterion  (e.g..  Ward  [13 ]>  Walker 
[l4].  Rice  and  Martin  [15]>  etc).  The  necessity  of  better  londerstanding  of 
these  ccmplex  problems  helped  motivate  the  present  study  of  the  loudness  of 
impulsive  sounds. 

2.2  Mechanism  of  Hearing 

The  perception  of  sovind  by  the  himnan  ear  is  a complicated  process  that 
is  understood  in  the  main,  but  has  aspects  that  remain  controversial. 

The  ear  consists  of  three  main  parts:  the  outer  ear  which  "matches" 
the  impedance  of  the  ear-drums  to  the  air,  the  miiddle  ear  which  transmits  the 
mechanical  vibrations  of  the  ear-drum  to  the  third  part  called  the  inner  ear. 
Inside  the  inner  ear  along  the  basilar  membrane  of  the  cochlea  the  sounds  are 
analysed  and  transferred  to  the  nervous  system. 

The  human  hearing  mechanism  acts  together  as  a microphone,  a highly 
selective  frequency  analyzer,  a sound  localizer,  and,  via  the  neurological  system 
and  brain,  an  interpreter  of  the  loudness,  the  pitch,  and  the  timbre  of  sound. 

The  acoustic  pressures  over  which  the  human  ear  can  operate  cover  the  range 
from  1000  microbars  down  to  0.0001  microbar,  which  is  extremely  wide.  The 
frequency  range  covers  about  ten  octaves,  or  approximately  from  20  to  20,000 
Hz  for  youthful,  healthy  ears. 

The  oldest  and  perhaps  most  comprehensive  theory  of  auditory  analysis 
was  develcped  by  Helmholtz  in  1863  [l6].  Since  then,  however,  the  Helmholtz 
assumptions  were  proved  doxzbtfvil  and  many  newer  theories  have  been  developed. 
Helmholtz  assumed  that  the  ear  separates  a complex  sound  into  sinusoidal 
components  corresponding  to  those  in  a Fovtrier  analysis  and  that  every  discrimi- 
nable  pitch  corresponds  to  one  particular  group  of  nerves.  These  groips  were 
svpposed  to  be  connected  to  the  specific  cochlear  segments  resonating  to  the 
specific  tones.  One  of  the  problems  not  satisfactorily  solved  by  this  theory 
refers  to  the  response  to  several  simultaneous  tones,  as  in  a chord:  this  is 
heard  as  a single  complex  tone,  and  is  not  decooposed  into  the  individual  pure 
tones.  It  is  also  hard  to  explain  the  behaviour  of  the  cochlear  "resonators" 
in  Helmholtz'  model. 

In  this  century,  von  Bekesy  (e.g.,  [17],  [I8],  [19])  stands  as  a giant 
in  auditory  research.  Using  surgical  and  physical  methods  on  test  animals,  he 
advanced  very  far  in  clarifying  the  role  of  the  tapered  helical  organ  called  the 
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cochlea.  By  means  of  a mechanical  model  he  gave  credibility  to  the  notion  of 
the  spatial  mapping  of  frequency  response  along  the  cochlea  ("place  theory"): 
the  lower  the  frequency,  the  closer  the  point  of  peak  vibratory  response  to  the 
apex  (small  end)  of  the  cochlea. 

In  general,  two 'frequency  theories" explaining  pitch  discrimination  were 
developed.  The  first  one  assumes  that  in  the  low  and  middle  ranges  of  the  toned 
spectrum,  pitch  is  determined  by  the  frequency  of  the  neural  in^piilses;  the  second 
one  assumes  existing  preliminary  mechanical  analysis  at  the  basilar  membrane. 

Both  of  the  classical  theories  (i.e.,  freqioency  and  place  theory)  presvq)pose 
that  the  frequency  analysis  is  oon5)leted  before  any  action  in  brain  can  take 
place,  and  both  have  many  theoretical  problems  to  account  for. 

Some  researchers,  like  Weddell  [20]  and  Pfafftaan  [21],  in  connection 
with  cutaneous  and  gustatory  senses,  produced  evidence  that  the  discrimination 
process  of  sensation  quality  depends  on  activity  in  a grovq)  of  nerves  rather 
than  in  one  nerve  fiber.  These  findings  s\;q)port  modern  views  on  auditory 
theory  suggesting  that  the  sensory  patterns  are  interpreted  by  the  bradn  and, 
therefore,  postulate  some  central  analysis  (Nordmark  [22]) . 

The  general  conclusion,  based  on  today's  knowledge,  is  that  there  is 
no  siniple  and  mique  relationship  between  physical  parameters  of  soiuid  and 
human  perception  of  sound.  From  the  subjective  point  of  view  the  physical 
characteristics  of  sound,  such  as  intensity,  frequency,  etc.,  maybe  described 
in  terms  of  three  basic  subjective  characteristics,  i.e.,  loudness,  pitch  and 
timbre.  None  of  them  is  dependent  strictly  iq)on  a single  physical  quantity, 
thovigh  it  is  possible  to  connect  each  of  them  primarily  to  some  physical  char- 
acteristic of  sound. 

The  sensation  of  pitch  is  determined  mainly  by  the  frequency  of  the 
soimd,  but  it  is  also  a function  of  the  intensity  and  wave-form.  Surprisingly, 
doiibling  the  frequency  does  not  double  the  pitch. 

Timbre  is  a subjective  characteristic,  which  makes  it  possible  to 
dlstii^gidsh  between  two  tones  having  different  wave-forms,  even  when  the 
intensity  level  and  fundamental  frequency  are  the  same;  but  it  is  also  a ftinc- 
tlon  of  intensity  and  frequency. 

Loudness  is  determined  by  the  intensity  of  the  sound  and  its  spectral 
distribution. 

2.3  Lovidness  of  Sound 

2.3.1  Definition  and  Msasurement  of  Loudness 

The  first  significant  atten^ts  to  define  and  measure  the  subjective 
quantity  of  so\md,  Associated  with  magnitude,  were  made  by  Fletcher  and 
Steinberg  [23]  in  1024  and  later  in  1925  by  Steinberg  [24].  In  1957  Fletcher 
and  Munson  [25]  defined  loudness  as  the  "magnitude"  of  sound  and  established  a 
1000  Hz  p\are  tone  as  a standard  tone  against  vdiich  the  loudness  of  other  sounds 
can  be  judged. 

Stevens  in  1936  [26]  introduced  the  widely  used  loudness  scale  with  a 
tmit  called  sone,  which  is  defined  as  the  loudness  of  a 1000  Hz  pure  tone  at  a 
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sound  pressxire  level  of  40  dB  (re  0.0002  micrbbar) . The  doubling  of  the  sound 
loudness  gives  a double  value  on  the  sone  scale.  It  was  mentioned  before  that 
loudness  as  a subjective  characteristic  of  souind  depends  iipon  the  intensity  of 
the  sound  and  also  is  a function  of  its  frequency  spectrum.  Various  investiga- 
tors determined  equal  loudness  contoxirs  for  pure  tones  and  for  band  noise. 
Fletcher  and  Munson  [25]  determined  the  sound  pressure  levels  of  pure  tones 
requilred  to  be  Judged  equally  as  loud  as  a 1000  Hz  reference  tone  over  nost  of 
the  auditory  frequency  range.  This  and  other  studies  are  not  entirely  in  agree- 
ment. A consensus  was  arrived  at  by  the  International  Organization  for  Standard- 
ization in  the  form  of  standard  ("nonnal")  equal-loudness  contours  for  pure 
tones  [27]-  The  ISO  recommendations  axe  commonly  used  as  the  statistical  result 
for  the  total  population  with  the  restrictions  for  the  position  of  a source, 
free  field  conditions  and  plane  wave  sound,  binaioral  listening,  etc. 

Recently  there  have  been  a few  attempts  to  extend  the  existing  equal- 
loudness contours  toward  the  lower  so-called  "infrasonic"  frequencies  (e.g.. 
Whittle  et  al  [28]). 

One  of  the  most  in5)ortant  problems  in  loudness  reseeirch  has  been  the 
development  of  a method  for  theoretical  prediction  of  its  level  for  the  complex 
sounds.  First  attempts  at  a loudness  calculation  based  on  the  frequency  spectral 
measurements  produced  very  conplicated  and  therefore  not  very  useful  procedures 
(Steinberg  [24],  Fletcher  and  Munson  [29]).  Two  of  the  newer  ones,  capable  of 
dealing  simply  with  broad  band  continuous  spectra,  were  adopted  by  the  Inter- 
national Standard  Organization  and  axe  widely  used  [30]. 

From  his  original  procedure  [31 ]»  [32]  Stevens  developed  a new  method 
called  the  Mark  VI  [33]  for  predicting  loudness  levels  of  noise  measured  in 
octave,  one-half  and  one-third  octave  bands.  This  procedure  adds  the  fraction 
of  so-called  "loudness  indexes"  in  the  bands  to  the  "loudness  index"  of  the 
loudest  band: 

Loudness  = S + f(ZS  - S ) 
m ' jn' 

where  S = sxaa  of  loudness  indices  of  all  bands, 

S^  = greatest  of  loudness  indices  in  any  one  band, 

f = fractional  portion  dependent  on  bandwidth. 

Values  axe  found  from  the  graph,  which  is  slightly  different  from  the  equal- 
loudness (sone)  contours  and  v&s  found  by  successive  approximation  (Stevens 
[33]) • The  resiiLt  of  the  calculations  is  obtained  in  sones;  however  the 
loudness  is  very  often  expressed  in  logarithmic  terms,  called  "phons" . The 
phon  is  eqviivalent  to  the  decibel,  but  instead  of  the  usual  definition,  similar 
to  the  decibel  definition  (20  logio  of  the  ratio  of  two  loudness  levels)  it  is 
derived  from  psychologicsil  units.  The  loudness  of  any  tone  on  the  same  equal- 
louidness  contour  has  a valxxe  in  phons  which  is  equal  to  the  sound  pressure 
level  in  dB  of  the  1000  Hz  tone  on  the  same  equal-lovidness  contour.  An  increase 
in  loudness  level  of  10  phons  is  approximately  equivalent  to  dovibling  the  subjec- 
tive loudness  (in  sones)  in  the  range  of  mid  and  high  levels. 

In  1971  Stevens  pviblished  a report  with  a new  modified  procedxjre 
called  Mark  VII  [34],  which  gives  the  perceived  level  of  loudness  or  noisiness 
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in  "PLdB"  units.  The  new  different  standard  reference  sound  is  used,  and 
the  procedure  allows  the  possibility  of  calculating  the  perceived  levels  of 
loudness  for  transient  sounds  (sonic-boom),  by  incorporating  the  Johnson- 
Robinson  procedvtre  [2]. 

The  ISO  also  endorsed  a procedure  developed  by  Zwicker  [35],  [36]. 
Together  with  Flottorp  and  Stevens  [37],  Zwicker  demonstrated  the  existence  of 
2k  so-called  "critical  bands"  in  the  audible  frequency  range.  Beyond  each 
critical  bandwidth  any  further  increase  in  a band  of  noise  was  found  to  have 
no  influence  on  the  amount  of  masking  produced  by  that  noise  on  a pure  tone  at 
the  centre  of  a band.  Zwicker  assumed  the  existence  of  a ftmctionsl  corre- 
spondence between  masking  and  loudness,  and  on  this  basis  developed  a method 
for  calculating  the  loudness  of  complex  sounds.  He  prepared  several  graphs 
on  which  the  vertical  divisions  are  marked  in  sones  and  the  horizontal  divisions 
divide  the  scale  according  to  the  equal  critical  bands  in  Hz  (approximated  by 
one-third  octave  bands) . Plotting  the  sound  spectrum  on  such  a prepared  graph 
and  drawing  in  the  lines  for  spread  of  masking  encloses  an  area  which  is  propor- 
tional to  the  total  loudness  of  the  sound.  Zwicker  defines  as  1 sone  the  area 
encompassed  on  the  graph  by  a one-third  octave  band  centred  at  1 kHz  at  a sound 
pressiire  level  of  40  dB  (including  the  correction  for  masking) . The  Zwicker 
and  Stevens  methods  give  estimates  of  loudness  levels  for  the  same  sounds, 
differing  by  approximately  3 to  5 phons. 

It  is  worth  mentioning, at  this  point,  the  existence  of  another 
fesponse  concept  called  perceived  noisiness.  The  perceived  noisiness  is 
defined  as  a subjective  iii5)ression  of  the  annoyance  or  unacceptability  of  a 
sound.  The  units  of  the  noisiness  are  called  "noy"  and  "PNdB"  in  analogy  to 
sone  and  phon. 

2.3.2  Loudness  of  Transient  Sounds 


In  a previous  peiragraph  we  discussed  the  general  concept  of  loudness 
in  the  context  of  steady- state  sounds.  In  the  case  of  transient  sounds  the 
situation  becomes  much  more  coii5)licated. 

Numerous  experiments  have  established  that  the  ear  responds  to  sound 
energy  as  integrated  over  a certain  time.  The  loudness  of  the  short  pulses 
grows  with  their  duration,  and  some  evidence  shows  that  the  loudness  level 
(which  is  logarithmic)  is  nearly  directly  proportional  to  the  logarithm  of 
the  duration,  for  sovinds  with  narrow-band  spectra.  In  general  we  can  regard 
the  whole  auditory  system  as  an  almost  linear  energy  integrator  for  transient 
sounds;  however  some  experiments  have  shown  different  behaviour  (stiller  [38], 
von  Port  [39])  8*  levels  near  the  threshold  of  hearing. 

An  early  assumption  concerning  temporal  integration  was  that  the 
ear  migpit  perform  as  a perfect  integrator  and  so  woiild  show  3 dB  reduction  of 
the  auditory  threshold  for  each  doubling  of  a signal  duration  (Gamer  and 
Miller  [4o]).  The  very  short  signals  as  well  as  very  long  ones  were  found  to 
deviate  from  this, simple  rule;  therefore  there  were  attempts  to  describe  the 
relationship  with  a three  line  graph  (Green  et  al  [4l]) . Harris  et  al  [42] 
in  1958  introduced  "critical  duration"  defined  as  the  intersection  of  the 
most  linear  portion  of  the  integration  function  with  the  abscissa  10  log(lt/l^ 
= 0 where  I^  = threshold  intensity  at  t = w. 


! 


I 


V 


In  1959  Pico?)  and  Bovcnan  £43]  derived  two  hypothetical  response 
functions  based  on  a siii?)le  electrical  model;  they  found  that  the  following 
equation  provides  a best  description  of  their  experimental  results  obtained 
near  the  threshold  of  hearing: 


10  iogjo(i/iJ  - -10  ioe^o(i 


(2.1) 


where  t = time  constant, 

t = duration  of  tone  pulse, 

I = threshold  intensity  of  tone  pulse, 

= threshold  intensity  at  t » op  (svistalned  tone) . 

Zwislockl  £44]  developed  a quantitative  psychophysiological  theory  which  gives 
an  explanation  of  the  phenomenon  of  the  louchiess  level's  dependence  upon  the 
stimulus  duration.  His  theory  is  based  on  the  evidence  that  the  appeirent 
teinporal  siamnatlon  of  acoustic  energy  is  a result  of  a neural  summation  at  a 
high  level  of  the  auditory  system. 

Due  to  the  large  differences  in  empirical  measurements  caused  by  the 
lack  of  accurate  e:q)erimental  conditions,  different  types  of  stimuli  used  in 
ejqierlments  and  large  individual  <31fferences  between  siibjects,  there  is  con- 
siderable disagreement  about  the  value  of  the  time* cons teuit  of  the  ear.  The 
reported  results  range  from  20  to  200  ms,  which  can  produce  up  to  10  dB  differ- 
ence in  loudness'.  The  first  investigations  of  this  problem  were  done  by  Bekesy 
£17]  in  1929 > who  obtained  the  value  of  the  constant  equal  to  I80  ms,  and  later 
by  Munson  £45]  with  the  resxilt  of  250  ms. 

Von  Port  £39] > using  gated  noises  of  various  bandwidths,  centre 
frequencies  and  levels  averaged  the  results  of  the  measurements  to  70  ms. 

Niese  £46],  sinanarlzlng  a nimber  of  studies,  obtained  the  average 
value  of  25  ms.  Other  e:q>eriments  ga^^e  (different  results:  Stevens  and  Hall 
£47]  about  150  ms,  Zwicker  £48],  £49]  100  ms,  [50]  from  10  to  20  ms, 

etc. 

In  all  cases  the  ranges  of  Judgements  between  individual  subjects 
were  reported  to  be  very  large,  whlcdi  probably  reflects  the  difficulties  of 
the  Judgements  themselves  or  substantial  individual  variance. 

Some  of  the  investigators  have  found  the  value  of  the  time-constant 
to  be  dependent  on  a sensation  level  (Small  et  al  £^0]) , while  others  did  not 
obtain  any  Indication  of  such  dependence  (e.g.,  von  Port  £39])*  The  controversy 
vfticther  the  time- constant  of  the  hvonan  au(U.tory  system  depend  vpon  the  frequency 
of  a stimulus  is  also  not  resolved.  The  loudness  of  transient  stimuli  can?>ared 
to  continuous  sowds  depends  also  upon  other  characteristic  parameters  describing 
the  waveform  and  its  duration.  The  rise  time  and  repetition  rate  are  strongly 
correlated  with  the  loudness  level  of  the  impulsive  sounds. 
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The  Influence  of  rise  time  on  the  loudness  has  been  studied  by  Cairter 
[51]  who  used  trlangulax  pulses  with  very  short  rise  times  (0.05  - 0.5  ms) . 
GJaevenes  [52]  used  tone  bursts  with  rectangular  envelope  and  rise  times  between 
30  and  9^0  ms.  From  both  experiments  authors  concluded  an  Increeise  in  loudness 
with  a decrease  in  rise  time.  Carter  suggested  that  the  relationship  between 
the  lovidness  level  eind  the  rise  time  is  of  a logarithmic  type.  Gustafsson  [53] 
used  in  his  experiments  pulses  from  a transient  sound  generator  having  rise 
times  from  0.3  to  10  ms.  Results  showed  the  connection  between  loudness  and 
rise  time  in  a form  of  a power  fiaiction,  which  tends  to  confirm  Carter's  results. 

The  effect  of  repetition  rate  on  the  loudness  of  short  pulses  was 
also  studied  by  Carter  [51J.  He  presented  data  showing  that  for  the  triangular 
transients  with  0.5  ms  rise  time  and  1 ms  duration,  having  the  repetition  rate 
from  1 pps  (pulse  per  second)  to  128  pps,  the  doubling  of  the  repetition  rate 
gives  3 dB  rise  in  loudness  of  the  signatures.  A repetitive  train  of  noise 
bursts  (say  10  ms  on  and  10  ms  off)  is  Judged  to  be  about  3-^  louder  than  a 
spectrally  similar  continuous  noise  of  equed  peak  level  (Pollack  [5^]);  however, 
the  continuous  sound  contains  about  8 dB  more  energy.  Zwicker  [48]  showed  that 
this  effect  can  be  replicated  when  subjects  axe  instructed  to  Judge  annoyance 
but  he  found  no  such  effect  when  they  were  instructed  to  Judge  loudness. 

2.4  Human  Response  to  Sonic  Boom;  Loudness 

2 .4 .1  Introduction 


Sonic-boom  research  involves  a whole  spectrum  of  different  fields 
and  problems.  The  research  in  physiology  covers  sleep  losses,  thresholds  of 
waking  and  stress- induced  responses  of  the  human  beings  and  animals.  From  a 
psychological  point  of  view  startle  measurements  and  task  interference  are  the 
main  subject  of  the  investigations.  Group  actions,  ccoplalnts  and  ccmmunlty 
protests  axe  paxt  of  the  sociological  programs. 

The  chosen  environment  and  investigative  method  in  sonic-boom  research,  , , 
as  well  as  in  research  involving  other  types  of  noise,  often  depends  upon  the  ^ < 

purpose  of  an  experiment  and  problems  connected  with  it.  In  laboratory  studies 
with  simTilated  or  reproduced  signals,  stimulus  parameters  and  the  investigative 
environment  situation  can  be  easily  controlled.  The  most  publicized  sonic-boom 
experiments,  however,  have  been  carried  out  \jnder  field  conditions,  using  actual 
aircraft  overflights.  Such  relatively  uncontrolled  field  studies  are  the  only 
way  to  do  community  research. 

2.4.2  Simulation  Techniques 

The  laboratory  research  on  sonic  boom  usvially  requires  long  and 
very  extensive  programs  with  a possibility  of  controlling  and  easily  varying 
the  sigpiature  parameters.  Therefore,  it  is  necessary  to  find  other  ways  of 
producing  test  signals  than  the  very  expensive  and  inconvenient  svpersonlc 
overflights.  With  this  motivation  veu’ious  types  of  simulators  have  been 
developed  in  a number  of  laboratories.  They  are  usually  highly  specialized  and 
based  on  different  principles  of  operation,  due  to  specific  research  requirements. 

One  scheme  of  sonic  boom  simulation  utilizes  a tiny  precision  model  of 
an  SST  placed  in  a wind  tunnel.  Projectiles  fired  in  a range  are  also  used. 

Propagation  experiments  may  be  conducted  with  these  simulators. 
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Shock-tube  or  con^ressed  air  driven  simulators  are  used  for  investi- 
. gating  propagation  phenomena  as  well  as  sonic  boom  response  of  hianans,  animals 
and  structures.  The  driver  is  coupled  to  a large  pyramidal  horn  down  which 
the  travelling  N-wave  propagates.  The  driver  device  may  be  a shock  tiabe  with 
diaphragm  (Warren  [55]}  Slutsky  and  Arnold  [56])  or  a specially  designed  com- 
pressed air  value  (Tamboulian  and  Peschke  [57]>  Glass  et  al  [^])«  A shock 
tiibe  with  a very  short  horn  formed  the  basis  of  a portable  sonic-boom  simulator 
for  field  studies  of  wildlife  (Gottlieb  [58]) . For  large-scalQ  outdoor  environ- 
mental experiments,  sonic-boom  waveforms  may  also  be  Simulated  by  multiple  or 
specially  shaped  (for  better  control)  explosive  charges. 

Another  type  of  simulation  facility  eii5)loys  a piston  as  driver.  In 
the  NASA  Langley  Research  Center  the  low-frequency  piston  facility,  with  elec- 
^ trdhydraullc  drive  and  cylindrical  test  chamber  were  eii?)loyed  for  structural 

and  environmental  testings  by  Edge  and  Mayes  [59]-  Similar  facilities  at  the 
Wright-Patterson  Air  Force  Base  and  at  the  Stanford  Research  Institute,  have 
been  used  specifically  for  testing  single  h\anan  subjects  (Lukas  and  Kr^er 
[60]). 

However,  the  most  convenient  and  flexible  technique  for  the  experi- 
mental testing  of  the  htonan  response  to  the  simulated  sounds  (e.g.,  sonic-boom) 
enploys  earphones  or  loudspeakers  for  the  simxalation.  In  order  to  obtain  the 
reqiolred  low  frequency  response  (down  to  DC  in  case  of  the  sonic-boom  spectrum), 
it  is  necessary  to  seal  either  the  earphone  or  the  loudspeaker- driven  chamber 
to  be  virtually  airtight.  Additionally  a specially  equalized  electronic  ampli- 
fier must  be  provided.  Sonic  boom  stibjective  experiments  with  the  application 
of  earphones  were  conducted  by  Ihtirner  {6l],  Ellis  [62],  Zepler  and  Harel  [1] 
and  later  by  Rood  [63]  at  the  University  of  Southampton.  This  type  of  facility 
reqioires  the  development  of  a device  which  will  secxue  the  airtight  fit  of 
earphones  to  the  head  of  a subject.  However,  the  main  disadvantage  of  the 
earphone  system  is  the  lack  of  the  ^diole  body  exposure  to  the  sonic  boom 
pressiire  field,  an  exposure  which  might  be  important  for  some  subjective 
studies. 

The  airtight  loiadspeaker-driven  chamber  overcomes  this  deficiency. 

In  1965  Pearson  and  Kryter  [64]  used  a 100  ft3  chamber  with  18  inch  loud- 
speakers mounted  in  a wall  and  driven  by  an  amplifier.  The  signals  were 
generated  by  means  of  the  device  giving  the  desired  electrical  waveforms  by 
following  a silhouette  placed  in  a "photoformer".  Also  in  1965,  a similar 
chamber  facility  was  developed  and  built  in  the  Bioacotistic  Laboratory  at 
Lockheed  for  research  by  Shepherd  and  Sutherland  [6].  The  loudspeakers  in 
this  facility  were  moianted  on  the  door  and  were  driven  by  high  power  DC  ampli» 
fiers  employing  special  equalizing  circuitry.  The  input  test  signals  were 
recorded  on  a magnetic  tape.  An  electronic  noise  squelch  circuit  was  employed 
' for  reducing  the  background  noise.  The  70  cubic  foot  chamber  allowed  over- 

pressures of  4.5  psf  for  long  rise  times  (10  ms),  and  2-3  psf  for  the  shortest 
rise  time  of  1 ms. 

j The  simulation  device  developed  in  the  University  of  Southampton 

' in  the  form  used  by  Rood  [63],  consists  of  the  120  cubic  foot  chamber  and  an 

N-wave  electronic  generator.  The  additional  frequency  compensating  network 
and  negative  feedback  system  allowed  the  attainment  of  1 psf  overpressure  at 
* 2 ms  rise  time. 
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The  current  investigation  was  carided  out  in  the  loudspeaker- driven 
simiiLation  booth  developed  at  the  Institute  for  Aerospace  Studies,  University 
of  Toronto,  which  is  described  in  Chapter  3,1.  Building  from  the  experiences 
with  the  earlier  booths,  siibstantial  advances  in  design  have  been  incorporated. 
Rise  times  as  low  as  0.2  ms  are  routinely  obtained  with  overpressures  of  order 
1 to  2 psf,  and  good  waveform  simulation. 

2.4.3  Prediction  Methods  for  Loudness  of  Sonic  Boom 

In  Chapter  2.3.1  we  discussed  methods  for  calctdating  the  loudness 
of  continuous  sounds.  Now  we  will  concentrate  on  the  methods  specifically 
developed  for  calcixlating  the  loudness  of  Inpulsive  sounds. 

In  1963  von  Port  [39]  suggested  a procedure  which  transformed  the 
energy  of  the  inpulsive  sotinds  to  such  a form  that  it  was  possible  to  apply 
the  equal-loudness  contours  for  steady-state  sounds  in  the  prediction  process. 
He  assumed  that  the  ear  integrates  the  energy  over  a 70  ms  interval.  The  value 
of  the  critical  time  was  found  through  subjective  measurements  of  the  loudness 
of  narrow  and  broad  bands  of  white  noise  pulses  as  a function  of  duration.  For 
pulses  shorter  than  70  m the  measured  sound  energy  was  divided  by  70  ms  and 
the  loudness  was  calculated  as  for  a continuous  sound.  For  repeated  ispulses 
von  Port's  procedure  is  more  cooplicated  and  involves  a different  decay  time 
(350  ms) . 


Niese  [46]  in  1975 » using  the  same  basic  ideas  as  von  Port,  proposed 
a method  which  differs  mainly  in  how  the  steady-state  measurements  should  be 
made.  The  method  starts  with  the  A- level  weigjited  calculation,  corrected  by 
the  size  and  frequency  of  the  largest  one- third  octave  band  measurement.  The 
integration  time  used  by  Niese  is  25  ms. 

Kryter  [65]  attenpted  to  develop  a mechanistic  calculation  procedure 
in  which  one  calculates  "perceived  noise  level"  (PNL)  over  0.5  ms  intervals. 
For  inpulsive  sounds  the  calculations  were  corrected  for  "startle". 


All  the  procedures  mentioned  above  were  designed  to  calculate  the 
loudness  of  the  inpulsive  sounds  in  general.  The  two  following  methods  were 
developed  with  the  particular  purpose  of  calculating  loudness  of  the  sonic- 
boom  signatures. 


Zepler  and  Harel  [1]  concluded  from  the  various  e:q)erimental  data 
that  the  mechanism  of  critical  bands  for  impulsive  sounds  is  not  clear.  There 
axe  some  indications  that  the  critical  bands  are  created  by  the  signals  them- 
selves or  that  the  parameters  of  the  mechanism  of  frequency  selectivity  vau-y 
under  intelligent  control.  Therefore,  they  iii5>lled  that  it  is  not  possible  to 
meastire  such  critical  bands  for  short  pulses.  In  their  research,  Zepler  and 
Harel  have  tried  to  find  a relationship  between  weighted  (in  relation  to  the 
equal-loudness  phon  contours) , energy  of  the  stimulJis  and  the  subjective 
loudness  of  a sonic-boom  slg;nature.  They  calculated  so-called  "weighted  energy 
density"  k|F(u))|  of  the  sound.  The  factor  "k",  which  converts  decibels  into 
phons,  is  ^osen  with  some  arbitrariness.  Based  on  the  enpirical  evidence, 
the  authors  eissumed  that  the  unit  step  fxinction  (at  127.6  dB  peak  pressure)  at 
the  frequency  of  50  Hz  affects  the  hearing  in  the  same  way  as  50  Hz  continuous 
tone  at  the  level  of  80  dB  SPL.  At  1 kHz  the  value  of  the  factor  "k"  is  unity. 
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The  general  relationship  is  as  follows: 


k - lo-Ao 


where  "n"  is  the  difference  between  phons  and  decibels  at  the  ^propriate 
level,  for  a particular  frequency.  In  their  report,  they  suggested  that  the 
reference  level  chosen  for  the  calculations  [1]  may  not  be  best  fitted  and 
Rood  [63]  has  con5)ared  his  experimental  resiilts  with  theoretical  calciaations 
based  on  this  method,  for  a few  reference  levels,  and  found  that  the  best  results 
were  obtained  for  a 100  phon  weighting  network.  The  chosen  level  detemines  the 
weighting  factor  "k"  The  subjective  loudness  is  proportional  to  the  area  under 
the  weighted  energy  density  ciarve. 

The  Zepler  and  Harel  method  for  predicting  loudness  was  also  svqiported 
by  experimental  evidence  by  Rice  and  Zepler  [66],  for  short  duration  waveforms 
caused  by  firing  a pistol  shot. 

The  second  sonic-boom  loudness  prediction  theory,  by  Johnson  and 
Robinson,  was  developed  in  two  papers,  published  in  1967  [67]  and  in  1969  [2]. 

It  follows  von  Port's  proposal  in  many  respects.  The  procedure  is  based  on 
the  assun^tlon  that  the  ear  Integrates  the  energy  of  a sound  arising  iii  a 
given  frequency  band  coiDpletely  if  the  whole  of  the  energy  is  received  within 
less  time  than  the  auditory  time  constant.  Thus,  the  loudness  of  the  short 
transients  is  determined  by  the  sound  energy,  not  by  the  sound  pressure  level , 
in  the  spectrum  bands.  The  value  of  70  ms,  as  a ccD^rcmise  between  the  growth 
and  decay  time  constants  of  the  integration  process,  was  chosen. 

The  energy  spectral  density  |f(u))  of  the  transient  stimuli  may  be 
obtained  through  digital  Fast  Foiarier  Transform  procedure  or  any  other  method 
suitable  for  the  short  sound  pulses.  The  Fovuier  transform  of  the  "ideal" 

N-wave  signature  determined  analytically  is  given  by  the  equation  (Ref.  [2]): 


T+t 
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where  T = diu*ation  of  the  signature, 

T = rise  time  of  the  signature. 


0)  (T  ^ t)T 


to  (T  ^ t)  t 


(2.2) 


In  order  to  sinplify  the  congputations , Johnson  and  Robinson  proposed  that  the 
energy  in  the  one-third  octave  bands  may  be  cstlculated  from  the  envelope  of 
the  spectrum,  since  for  the  frequency  range  of  interest  the  actual  spectrum 
oscillates  rapidly  (rate  1/ duration)  within  the  one- third  octave  bands  (see 
Fig.  1) . The  envelope  of  the  spectrum  function  is  defined  by: 
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because  the  mnyi Tna.  of  the  function  from  the  equation  (2.2)  occur  at  the  frequen- 
cies 0)  = 27m/T.  The  area  vinder  the  envelope  represents  approximately  twice  the 
reeil  band  energy,  therefore  the  energy  in  each  one-third  octave  band  is  propor- 
tional to; 

"2  sln^  ^ 

du)  (2.4) 

(>) 

J, 

for  an  N-wave. 

In  the  atten^jt  to  relate  the  impulse  energy  spectral  density  to  the 
steady-state  sound  loudness  calculation,  Johnson  and  Robinson  defined  the 
"effective  band  pressxire  levels"  of  sonic  boom,  as  to  be  equal  to  the  levels 
of  a continuous  steady  noise  which  over  a period  equal  to  the  critical  time 
contain  the  same  energy.  This  is  obtained  by  dividing  the  energy  spectrum 
density  levels  (calculated  directly  from  the  waveform  by  the  Foizrier  transform, 
or  produced  in  analog  form  by  the  spectrum  analyser)  of  the  transient  sound  by 
the  assumed  critical  time.  The  method  is  valid  only  for  pulses  shorter  than 
the  critical  time.  The  resulting  quantity  has  the  dimension  of  the  sound  power 
and  can  be  expressed  as  a sound  pressure  level  in  decibels,  std table  for  loud- 
ness calculation  by  any  of  the  existing  procedures  for  continuous  soiands  (see 
Chapter  2.3) . 

The  condition  for  the  pulse  duration  refers  to  the  effective  duration 
of  the  response  to  the  stimulus  at  the  level  of  the  cochlear  frequency- selective 
mechanism  and  not  to  the  physical  duration  of  an  ilIp^^lse.  Johnson  and  Robinson 
have  estimated  [67]  that,  for  typical  sonic  bangs  with  dviration  from  100  to 
250  ms,  even  if  the  response  for  the  lowest  frequency  bands  is  longer  than  70 
ms,  the  error  will  be  less  tlian  1 dB.  This  is  based  on  the  assxmption  that  the 
human  frequency-analytic  mechanism  has  the  properties  approximated  by  a set  of 
the  one-third  octave  filters  at  the  higher  frequencies  and  similar  set  of  filters 
with  wider  bandwidth  than  one-third  octave  at  the  lower  end  of  the  audible  spec- 
trum. It  can  be  assumed  that  for  these  filters  the  significant  part  of  the 
response  for  the  brief  pulse  is  of  the  order  l/Af  (Af -bandwidth) , which  means 
that  the  sonic-boom  N-wave  signatures  are  perceived  by  the  auditory  system  as 
two  separate,  individual  ptilses,  roughly  less  than  100  ms  long  for  all  bands 
below  the  50  Hz  band,  approximately. 

Having  the  one-third  octave  band  pressure  levels  in  decibels 
(re  2 X 10“5  N/itf^) , a few  corrections  are  introduced  before  applying  the  loud- 
ness calculation  procedure.  First,  because  the  duration  of  the  typical 
sonic-boom  signature  is  long  enough  so  that  the  two  bangs  (front  and  rear) 
are  heard  separately,  the  calculated  band  pressure  levels  are  reduced  by  3 dB . 

It  is  not  clear,  however,  what  the  limitations  for  this  correction  are  in  the 
sense  of  the  range  of  durations  for  which  it  should  be  applied.  Secondly,  in 
order  to  take  into  account  the  high  energy  levels  at  very  low  frequency  (e.g., 

3-4  Hz) , each  beind  pressure  level  below  the  50  Hz  one-third  octave  band  is 
reduced  by  means  of  the  equal-loudness  contours  to  that  level  at  50  Hz  which 
produces  the  same  loudness  level  in  phons.  These  weighted  levels  are  then 
combined  with  the  existing  level  in  the  50  Hz  band  by  simple  decibel  addition. 

The  data  prepared  this  way  may  be  used  for  loudness  evaluation  based  on  eoiy 
standard  procedure.  Johnson  and  Robinson  have  recommended  the  Stevens'  Mark 
VI  procedure  [33] > because  of  its  great  simplicity  and  the  good  results 
achieved  by  its  \ise. 

Although  Illustrated  in  the  foregoing  for  sonic-boom  N-waves,  the 
Johnson-Roblnson  predictive  method  has  a much  broader  affiliation:  it  may  be 
used  to  calcvilate  the  loudness  of  arbitrary  Inpulslve  wavefonos.  The  step-by- 
step  handbook  procedure  for  our  conputer.  implementation  is  detailed  in  Appendix  8. 
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2.4.4  Investigations  of  Sonic  Boom  Loudness:  Review 


A few  years  £igo  research  into  the  acceptability  of  the  sonic-boom 
by  hviinan  beings  was  mainly  concerned  with  possible  tengjorary  threshold  shift 
(TTS)  or  permanent  threshold  shift  (PTS)  in  hearing  response.  It  was  con- 
cluded, however,  that  these  were  ne^igible,  so  that  the  sonic-boom  can  be 
disregarded  as  a threat  to  the  auditory  system  (e.g.,  von  Gierke  and  Nixon 
[68])  Today  the  primary  direction  of  the  research  is  concentrated  on  the 
iii?)ressions  of  loudness  and  annoyance. 

Field  studies  of  sonic-boom  have  been  conducted  extensively  in  the 
USA,  England,  France,  and  Sweden.  Coii5)arative  investigations  of  loudness  and 
annoyance  for  sonic-boom  versus  jet  aircraft  flyover  noise  were  carried  out 
in  Great  Britedn  in  Project  Westminster  [69]  and  in  the  USA  in  the  Edwards 
AFB  Program  [70].  In  both  cases  the  observations  were  made  indoors  and  out- 
doors. A representative  resvilt  may  be  quoted  from  the  flight  tests  at  Edwards 
AFB , " . . . Forty  percent  of  those  polled  rated  booms  of  about  the  level 
projected  for  current  svq)ersonic  transports  in  cniise  (30  psf)  as  unacceptable 
for  outdoor  listening;  for  indoor  listening  the  corresponding  percentage  was 
30^(."  [71]. 

Johnson  and  Robinson  [67]  in  1967  conducted  a series  of  field  experi- 
ments in  which  6l  subjects  made  judgements  of  the  relative  annoyance  of  sonic- 
boom,  explosions  and  jet  aircraft  noise,  using  the  method  of  direct  magnitude 
estimation.  The  resvilts  of  the  subjective  experiments  correlated  closely  with 
calculated  loudnesa  levels,  suggesting  the  possibility  of  application  of  the 
proposed  theory  for  different  types  of  spectra.  The  correlation  with  perceived 
noise  and  peak  overpressure  was  poor.  In  their  later  paper  [2]  the  prediction 
theory  for  sonic-boom  was  used  in  theoretical  calculations  of  loudness  level 
as  a fxmction  of  rise  time,  duration,  and  delay  time  (between  the  incident  and 
the  ground-reflected  waves) . The  authors  assumed,  taking  into  account  previous 
work,  that  the  difference  between  the  spectrum  of  the  real  recorded  sonic-boom 
and  an  Idealized  approximation  consisting  of  to  seven  straight  line  segments 
are  of  little  consequence  in  detennining  loudness.  Therefore,  the  calciU-ations 
were  made  for  the  analytically  derived  Fourier  transforms  of  the  idealized 
N-wave  signatures,  using  the  proposed  prediction  procedure.  The  results  show 
no  variation  in  the  loudness  within  0.1  phon  in  the  duration  range  from  100  to 
500  ms  for  an  N-wave  having  0.1  ms  rise  time  and  zero  delay  time.  The  change 
in  rise  time  at  grotuid  level  (i.e.,  delay  time  equal  to  zero)  from  16  to  0 ms 
gives  the  increase  in  loudness  of  25  phons.  Within  the  range  from  16  down  to 
4 ms  loudness  rises  at  a rate  of  about  3/4  phon  per  millisecond  increase  of 
rise  time,  and  below  4 ms  the  increase  is  even  sharper. 

The  influence  of  delay  time  (between  the  direct  and  ground- reflected 
waves)  depends  upon  the  rise  time  value,  but  in  general  the  loudness  declines 
with  increase  of  delay  time. 

Turning  now  from  field  studies  to  laboratory  studies,  at  the  Univer- 
sity of  Southajqpton  research  into  the  subjective  effects  of  the  sonic-boom 
has  been  carried  out  by  Ellis  [62],  Zepler  and  Heu'el  [1]  and  Little  [72]  using 
earphones,  and  by  Rood  [63]  vising  a loudspeaker- driven  booth. 

Zepler  and  Heurel  in  their  subjective  tests  used  specially  designed 
e£UT)hones  developed  previously  by  Thumer  [61].  The  frequency  response  of 


13 


J 


earphones  was  contrived  to  be  practically  flat  between  zero  and  1.5  kHz. 

This  flat  frequency  response  allowed  reproduction  of  N-wave  signatures  with 
relatively  good  quality,  but  the  performance  depended  strongly  upon  air 
leakage  which  was  hard  to  avoid.  The  maximum  pressure  obtainable  was 
approximately  2 Ib/ft^. 

The  comparison  technique  was  used,  with  a continuous  tone  of  UOO  Hz 
as  a reference  signal.  While  this  reference  stimulus  has  the  advantage  of 

giving  results  directly  related  to  the  definition  of  loudness  xinits,  it  was 

reported  to  be  very  difficult  for  the  test  observers  to  conpare  loudness  of 
such  extremely  different  signals.  Altogether  ninety  subjects  took,  part  in 

the  experiment  and  the  stsindard  deviation  for  the  last  series  of  the  tests  for 

20  svibjects  was  reported  by  Zepler  and  Harel  as  being  about  4.5  dB.  Despite 
the  large  discrepancy  between  subjects,  resialts  for  individual  subjects  were 
usually  consistent  within  a few  decibels. 

They  found  a decrease  of  the  loudness  level,  at  the  same  overpressvire, 
equal  to  l6  dB  within  the  rise  time  range  from  0.3  to  10  ms.  This  result  was 
confirmed  by  theoretical  calcvilations  based  on  the  new  proposed  method.  Both 
the  experimental  and  theoretical  estimations  have  shown  that,  for  rise  times 
larger  thsui  1 ms,  variations  in  rise  time  at  constant  overpressure  had  more 
effect  than  the  corresponding  variation  in  the  overpressure  at  constant  rise 
time. 

The  authors  performed  an  additional  test  with  high  pass  filtered 
N-wave  signatures  showing  no  loudness  effect  for  cut-off  frequencies  lower 
than  4o  Hz,  for  N-waves  having  rise  times  of  1 ms  and  3 They  also 
reported  that  when  the  duration  of  the  N-wave  signature  is  shortened  to  50  ms 
or  less,  the  two  bangs  heard  for  longer  signatvires  separately  merge  into  one 
one  with  an  increase  of  about  3 dB  in  loudness. 

In  1973  Rood  [63]  piiblished  a Ph.D.  thesis  in  which  he  extended 
Zepler  and  Harel 's  research.  He  measured  the  absolute  subjective  loudness  of 
half  N-wave  signals,  using  a loudspeaker- driven  booth  as  a simulation  facility. 
The  pressure  booth  enabled  him  to  produce  signals  with  rise  times  from  l6  down 
to  2 ms  and  overpressures  from  0.5  to  3-0  psf  (the  latter  for  the  longer  rise 
times) . The  experimental  techniqvie  involved  the  conparison  of  N-waves  with 
200  Hz  and  400  Hz  tone  bursts,  20  ms  long.  In  his  paper  Rood  reported  some 
difficulties  in  finding  the  effective  pitch  and  siibj active  duration  of  the 
N-wave  signatxire  to  be  matched  by  reference  signals . A conparison  between 
tests  with  earphones  and  the  pressure  booth  was  also  conducted,  but  no  signi- 
ficant difference  in  respect  to  subjective  loudness  determination  was  foxuid. 

Rood's  experiments  were  carried  out  for  three  different  overpressure 
levels:  0.5  psf,  1 psf  and  2 psf.  The  loudness  conparison  test  results  were 
converted  into  the  absolute  loudness  level  scale  thro\;igh  the  experimentally 
established  relationship  between  the  reference  tone  bursts  and  500  ms  long 
continuous  sotind  and  these  values  were  coopared  then  with  the  theoretical 
, predictions.  The  standard  deviation  varied  from  2.44  to  7.l6  dB.  The  best 

I fit  to  these  results  was  obtained  by  the  Zepler  prediction  procedure,  for  100 

phon  weighting  curve.  In  general,  the  change  from  2 ms  to  l6  ms  in  rise  time 
caused  a decrease  of  subjective  lotidness  level  of  17  to  19  phon,  depending 
( ipon  the  peak  pressure.  It  was  found  that  while  the  doubling  of  the  over- 

[ pressxire  results  in  an  increase  of  loudness  by  6 phons  for  longer  rise  times, 

it  is  marginally  less  for  shorter  rise  times. 
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On  CQn5)aring  his  results  with  earlier  works  hy  Zepler  and  Harel  [1] 
and  Ellis  [62],  Rood  concluded  that  the  previous  results  underestimated  the 
true  loudness  levels  by  12  phons  because  of  the  test  method,  which  involved 
con?)arisons  of  the  iii5)Ulsive  sounds  with  the  continuous  sound  reference. 

The  most  can5)lete  and  significant  experimental  work  was  done  in  1968 
by  Shepherd  and  Sutherland  [6]  at  the  Stanford  Research  Institute.  They  employed 
the  Lockheed  loudspeaker-driven  sonic-boom  simulation  booth  (see  Chapter  2.4.2). 
The  total  subjective  test  series  was  divided  into  three  parts;  in  each  the  paired 
con^jarison  teclinique  was  used.  One  boom  in  each  conparison  pair  was  designated 
the  standard,  the  other  being  the  boom  under  test;  the  relative  aaplitudes 
were  adjusted  for  equal  subjective  loudness.  The  relationship  between  relative 
anplitude  versus  rise  time  and  duration  was  defined  at  three  different  anplitude 
levels.  Values  of  the  parameters  of  sonic-boom  signatures  used  in  the  experi- 
ments ranged  from  1 to  10  ms  for  rise  time,  from  100  to  500  ms  for  duration 
and  from  0.8  to  2.4  psf  for  overpressure.  Similar  se^s  of  conparison  tests 
were  repeated  for  annoyance  as  a judgement  criterion. 

In  the  second  part  of  the  experiment  the  subjects  evaluated  the 
loudness  of  half  N-wave  signatures  with  additioneil  sawtooth  serration,;  Jn 
the  last  part  there  was  an  additional  sawtooth  having  different  inteipeak 
spacing.  Shepherd  and  Sutherland's  basic  results  indicate  no  significant 
change  in  loudness  (inferred  from  amplitude  ratio  in  dB)  of  the  half  N-wave 
with  respect  to  duration  while  the  increase  of  the  rise  time  from  1 ms  to 
10  ms  was  followed  by  a decrease  in  loudness  of  13  dB. 

Increasing  the  amplitudes  of  the  two  signals  of  a pair  in  the  same  . 
proportions  had  no  effect  on  the  relative  loudness.  Also  the  difference 
between  loudness  and  annoyance  judgements  were  insignificant,  except  for  those 
involving  signatures  with  10  ms  rise  time. 

The  addition  of  a spike  to  the  half  N-wave  to  bring  the  amplitude 
to  3 '3  psf  was  tried;  it  provided  around  8 dB  difference  in  loudness,  when 
compared  with  the  half  N-wave  signature  without  spike,  of  overpressure  1.6  psf. 

There  was  no  significant  difference  in  loudness  between  signatures 
with  added  sawtooth,  so  long  as  this  additional  fine  structure  did  not  change 
the  general  shape  of  the  waveform  and  modify  its  rise  time . 
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3.  EQUIPMEUI 

3.1  UTIAS  Sonic  Boom  Simulator  - General 

In  order  to  obtain  more  significant  and  precise  data  on  human 
response  to  sonic  boom  and  other  types  of  noise,  the  development  of  a loud- 
speaker-driven booth  was  initiated  in  1970  at  OTIAS.  A series  of  modifications 
and  in^jrovements  led  to  a capability  of  faithful  simulation  of  the  sonic  boom 
N-wave  and  other  transient  waveforms.  A paxticuleir  feature  was  the  attainment 
of  rise  times  as  short  as  0.2  ms  at  an  overpressure  of  1 psf  (longer  rise  times 
at  higher  overpressures  up  to  4 psf) . This  rise  time  is  five-  to  ten-fold 
shorter  than  that  of  the  earlier  cited  facilities  (Shepherd  ajid  Sutherland 
[6],  Rood  [63]). 

The  UCIAS  Sonic  Boom  Simulation  Booth  (Fig.  2)  consists  of  an  air- 
tight 2.1  m3  volume  chamber  driven  by  12  loudspeakers  in  two  sizes  mounted  in 
the  aperture  of  the  wall  faced  by  subjects.  The  booth  features  a double-wall 
plywood  construction  with  Inside  wall  surfaces  heavily  lined  with  soiaid 
absorbing  semirigid  fiberglass  material  to  minimize  high  frequency  reflections 
and  consequent  resonances;  the  free-air  volume  is  thus  reduced  to  about  1.3 
The  absorbing  material  was  chosen  by  taking  into  account  easiness  of  attachment, 
porosity,  absorber  thickness  and  absorption  coefficient. 

An  Eiirtight  window  was  installed  on  the  door  to  allow  visual  observa- 
tion of  the  subject  seated  inside.  In  the  original  form  the  ventilation  system 
was  built  in,  which  allowed  a very  slow  flow  of  fresh  air  into  the  booth  during 
operation.  It  was  found  impractical,  however,  and  was  later  removed  to  minimize 
air  leakage.  The  booth  is  also  equipped  with  an  internal  door  release  mechanism 
and  locking  system,  intercom  and  a support  system  for  the  microphone. 

3 -2  UglAS  Sonic  Boom  Simulator  Electronic  System 

The  UEIAS  Sonic  Boom  Simulator  ertploys  six  15  inch  low-frequency 
loudspeakers  (Altec  Lansing  woofers,  model  515B)  and  six  8 inch  medium- frequency 
loudspeakers  (Radio  Shack  30  watts,  type  40-1286  with  the  frequency  range  up  to 
20  kHz)  to  cover  the  total  desired  frequency  range  from  about  0.1  Hz  to  5 kHz.  The 
two  types  are  displaced  in  their  mounting  in  such  a fashion  as  to  equalize  the 
effective  travel  time  to  the  siabject's  head.  Thus  the  signals  eirrive  in  step, 
so  that  the  wave  front  of  the  conposite  is  virtually  as  sharp  as  that  from  a 
single  medium- frequency  loudspeaker. 

The  electronic  system  consists  of  a crossover  circuit,  an  equalizing 
network  and  fo\ir  100  W power  amplifiers.  The  block  diagram  of  the  total  system 
is  shown  in  Fig  3.  The  low  and  high  frequency  part  of  the  input  signal  are 
separated  by  the  crossover  circtolt-  The  (adjustable  equalizing  network  conpen- 
sates  for  the  major  part  of  the  speaker  and  booth  coloration  of  the  frequency 
response  and  hence  elixninates  much  waveform  distortion.  (The  block  diagram  of 
the  equalizing  circuit  is  shown  in  Fig.  4.)  The  main  part  of  the  equalizing 
network  is  a two  channel  Altec  Lansing  "Acousta-Voicette"  equalizer.  Model  729A, 
contedning  twenty-four  one-third  octave  filters  centred  at  frequencies  from 
12,500  Hz  to  63  Hz.  Each  filter  is  adjustable  over  a range  ± 12  dB.  To  carry 
on  below  63  Hz  it  was  necessary  to  include  a series  of  low  pass  filters  in 
parallel  in  the  low  frequency  channel  (the  Altec  unit  cuts  off  below  8 Hz). 

The  total  signal  in  this  channel  is  summed  in  a specisil  circuitry. 

Carefiil  adjustment  of  these  filters  in  both  high  and  low  frequency 
channels  made  it  possible  to  obtain  a relatively  smooth  frequency  response  of 
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the  total  system.  This  is  supplemented  by  a special  noise- squelch  circuit 
consisting  of  an  offset  eliminator  and  a noise  eliminator  system  to  decrease 
the  background  noise.  The  system  is  triggered,  permitting  the  electronic 
network  to  be  switched  off  during  the  silent  intervals  between  the  test  pairs. 

The  output  signals  from  the  equalizing  network  axe  an^Jlified  by  four 
100  W power  amplifiers  with  nearly  flat  frequency  response  down  to  DC  and  less 
than  3 dB  drop  at  100  kHz . 

3.3  Generation  of  the  Input  Signals 

As  the  original  input  source  for  tbi  UPIAS  Sonic  Boom  Simulator,  an 
elaborate  electronic  function  generator  was  developed.  This  analog  facility 
had  the  capability  of  fitting,  within  limits,  an  arbitrary  voltage-time  signa- 
ture by  means  of  100  straight-line  segments  with  slopes  adjustable  by  individual 
knobs.  The  a4justme.it  turned  out  to  be  very  laborious;  but  worse  than  that 
the  signatures  were  not  stable  in  time  because  of  a problem  with  voltage  drift 
and  residual  noise  (Gottlieb  [33]) • Thus  the  function  generator  was  abandoned 
in  favour  of  a coii5)Uter-based  approach. 

The  test  signals  used  in  the  current  project  were  generated  from 
mathematical  formulas  by  the  HP  2100A  digital  con^juter  with  24k  memory,  magnetic 
tape  unit  and  conventional  D/A  and  a/D  equipment.  The  basic  block  diagram  of 
the  generating  sequence  is  shown  in  Fig.  5* 

A sin^jle  FORTRAN  program  generated  the  series  of  test  pairs  each 
having  a different,  randomly  presented,  overpressure  ratio  - according  to  the 
supplied  data.  The  whole  series  was  recorded  on  magnetic  tape  in  digital  form, 
and  then  after  D/A  conversion  recorded  in  analog  form  (Bruel  and  Kjaer  two 
channel  FM  tape  recorder  type  7001) . The  program  working  with  the  D/A  converter 
also  generated  the  triggering  signal  recorded  on  the  second  channel  of  the  tape 
recorder . 


As  mentioned  before,  the  equalizing  network  can,  within  its  limitations, 
reshape  the  frequency  response  of  the  simulation  system.  However,  the  adjustment 
is  not  sufficient  to  completely  eliminate  distortions  introduced  by  the  system 
and  achieve  a flat  frequency  response.  In  a case  of  more  conplex  waveforms 
(e.g.,  minimized  sonic-boom,  "flat-top"  signatures),  the  distortions  due  to  the 
non-ideal  frequenoy  characteristic  of  the  system  become  more  critical,  signifi- 
cantly changing  the  general  shape  of  a waveform.  Therefore,  a novel  ccmputer- 
aided  method  for  countering  the  distortion  was  developed.  Using  the  measured 
transfer  function  of  the  whole  simulation  system  (i.e.,  electronic  circuitry, 
loudspeakers,  booth)  the  input  signatures  axe  "predistorted"  in  frequency  domain 
to  counter  the  simulator  distortion.  The  theory  behind  this  counter- distort! on 
approach  is  given  below. 

Assuming  that  our  simulation  system  is  linear  and  time  invariant,  the 
output  y(t)  for  any  arbitrary  input  x(t)  is  given  by  the  equation: 

00 

y(t)  =J  h(T)  x(t  - T)dT  (3.1) 

-00 
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For  a physically  realizable  system  h(T)  = o for  t < 0,  where  h(T)  is  defined 
as  a system  response  to  the  unit  iii?)ulse,  applied  a time  t before. 


The  Laplace  transform  of  h(T)  gives  a transfer  function  of  the 
constant  parameter  linear  system  H(s),  defined  as: 


H(s) 


(3.2) 


where  s = 7 + Ju, 


The  dynamic  characteristics  of  such  a physically  realizable  and 
stable  system  can  be  described  by  a frequency  response  fvaiction  H((i))  which  is 
defined  as  a Fourier  transform  of  h(  t)  : 

00 

H(U))  = Jh(T)  dT  (3.3) 

o 

In  a special  case,  when  7=0  both  expressions  (i.e.,  3*2  and  3*3) 
are  equivalent  auid  the  frequency  response  function  may  replace  the  transfer 
function  with  no  loss  of  useful  information  [73].  Applying  the  Fovurier  trans- 
form to  equation  (3-1)  we  get  as  a result: 


Y(w)  = H(w)  X(u) 


(3.4) 


where  H(u))  = |H(w)  |e  ^ is  a transfer  fimction  of  the  system  and  Y(u)  and 
X(u)  are  Fovurier  transforms  of  the  output  auid  input  signals  respectively.  In 
the  ideal  case,  when  the  system  is  linear  and  its  parameters  ea*e  ccaastant  in 
time,  the  transfer  function  K(u))  is  a function  of  a frequency  only.  The 
xisefulness  of  equation  (3-4)  lies  in  the  fact  that  the  transfer  fxmction 
concept  allows  one  to  separate  the  system  characteristic  from  the  input  signal. 

We  may  now  represent  the  Sonic  Boom  Simulator  in  the  form  of  a black 
box  with  the  conplex  transfer  function  r(u).  If  we  put  an  ideal  N-wave  signature 
into  the  system,  we  would  come  ip  with  the  oirtput  signal  given  by  the  following 
eqiiation  (see  the  outline  of  the  predistortion  scheme  in  Fig.  6): 

F2(w)  = F^(u)  r(u>)  (3.5) 

where  F^(u)  = Fourier  transform  of  the  input  signal  (e.g.,  N-wave  signature), 
F^iu)  = Fourier  transform  of  the  resxU.ting  distorted  output  signal. 

Knowing  the  Fourier  transform  of  the  input  signal  we  can  recover  the 
transfer  function  of  the  simulating  system  from  the  spectrum  of  the  output 
signal: 

Km 

r(«)  = r--  (3.6) 

F^(u) 


f 
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The  inportant  thing  in  this  approach  is  that  it  must  be  insured  that  an  input 
signal,  measxiring  the  transfer  function,  will  excite  the  system  over  the  entire 
frequency  range  of  interest. 

Dividing  the  ccn^ilex  spectrum  function  of  the  desired  waveform  f^(w) 
by  the  transfer  function  of  the  simulating  system,  given  by  equation  (3.6),  will 
yield  a new  "predistorted"  Fourier  transform  function: 


^4(“)  = ="  F (u.)  

^ r(u.)  3 


(3.7) 


After  transforming  back  to  the  time  ’dcanain,  the  new  waveform  Fj.(t) 
fed  back  into  the  system  will  counter  the  distortions  of  the  non-ideal  transfer 
function  producing  the  desired  output  signal  F3(t) . In  terms  of  the  Fourier 
transform  functions  the  total  scheme  can  be  described  by  the  following: 


= -2 — TV 
r(o,) 


Output,  F^(u)  = Fjj((.))  r(u)  = ^ r(u)  = F^Cu),  desired  output 


(3.8) 


The  above  procedure  has  been  computerized;  the  entire  process  is  done 
by  a single  FORTRAN  program  (see  Appendix  7)  in  the  HP  2100A  computer  xising 
Fast  Fourier  Tretnsform  hardware.  The  computer  generated  "predlstorted"  signa- 
tures are  recorded  on  the  Bruel  and  Kjaer  FM  tape  recorder  and  played  back  into 
the  amplifiers  of  .the  Sonic  Boom  Simulator  in  the  sane  fashion  as*for  simply 
generated  signals.  The  simplified  block  diagram  of  the  generating  urocess  is 
shown  in  Fig.  7. 


From  the  theoretical  point  of  view  the  predistortion  scheme  should 
give  as  a result  ideally  clean  pressure  signatures,  identical  to  the  desired 
waveforms,  inside  the  booth.  This  is  never  achieved,  however,  in  practice. 
Beginning  with  the  simmlator  system,  we  made  the  assumption  of  linearity,  which 
is  only  an  approximation  of  reality.  Mainly  because  of  the  speakers  the  simula- 
tion system  exhibits  some  nonlinear  effects;  these  limit  the  accuracy  of  the 
predistortion  procedure  because  the  effective  transfer  function  of  the  system 
(?*)  is  also  a function  of  the  input  signal  itself.  The  transfer  function  depends 
also  upon  timie,  which  is  ca\xsed  by  changing  paramaeters  with  time  (e.g.,  air 
leakage) . 


The  discrete  Fast  Fourier  Transform  used  in  the  predistortion  proce- 
dure is  in  a sense  inexact;  it  is  only  an  approxinatlon  to  the  continizoi^s 
Fourier  transform  by  virtue  of  errors  due  to  discrete  sampling  in  time  and 
trvincatlon  of  the  sampled  signals.  Each  of  these  effects  can  be  miinlimized, 
but  the  lirndted  memory  of  the  ccaputer  and  coeputing  time  limitation  forces  a 
programmer  to  choose  a suitable  cooprcnlse  (more  details  about  FF7  are  given 
in  Chapter  3.4.2) . 

The  specific  characteristics  of  the  spectr;sm  fvinctions  of  the  signals 
used  in  this  investigation  caused  further  problems  for  the  predistortion  proce- 
dure. The  procedure,  it  will  be  rec6d.7.ed,  requires  the  division  of  one  spectnmm 
by  another.  However,  the  spectra  of  li^ulsive  signatures  have  midtiple  zeros 
(e.g.,  the  zeros  of  an  H-wave  spectnim  are  spaced  at  intervals  of  1/ duration); 
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thus  the  division  leads  to  O/O  at  these  points.  A number  of  different  schemes 
for  avoiding  divisions  by  zero  were  tested;  the  one  chosen  in  the  final  program 
substitutes  the  values  for  the  ideal  transfer  function  (i.e.,  the  constant 
g£iin  transfer  function  which  does  not  distort  the  input  aignal)  at  points 
where  the  real  transfer  function  cannot  be  defined  because  of  the  division-by- 
zero  situation  (see  equation  3-6). 

This  arbitrary  scheme  for  avoiding  the  divisions  by  zero  may,  of 
course,  yield  errors,  but  the  problem  was  fovind  to  be  non-critical  because  of 
the  limited  nunber  of  these  zero  points  in  practice  (i.e.,  the  discrete  points 
calculated  by  FFT  procedure  may  not  match  the  frequencies  of  zeros) . An 
additional  means  of  reducing  the  effect  of  O/O  was  afforded  by  using  the  same 
signals  for  testing  the  transfer  function  as  those  to  be  predistorted,  because 
zero  points  have  the  same  position  in  both  spectra.  It  eilso  assures  the  proper 
freqviency  range  of  the  testing  spectrum  and  reduces  effects  of  nonlinearity  of 
the  system. 

All  the  problems  discussed  above  introduce  errors  in  the  con5)uted 
spectrum  of  the  predistorted  input  signal;  these  are  reflected  in  the  time 
domain,  often  in  the  form  of  "ringing".  In  Fig.  8 there  are  shown  some 
examples  of  the  input  electrical  "predistorted"  signals  (top)  and  the 
resulting  output  pressure  signatures  recorded  inside  the  booth  (bottom) . 

Figttre  9 shows  the  comparison  between  pressure  signatures  recorded  inside  the 
booth,  without  and  with  predistortion.  The  latter  figure  illustrates  the  sig- 
nificant iii5)rovement  in  the  quality  of  the  waveforms  for  the  predistorted 
signals.  The  remaining  fine  structure  is  due  to  the  procedure  errors. 

Figures  10  and  11  show  the  respective  energy  spectrum  density  func- 
tions obtained  through  the  Fast  Fourier  transform  procedure  for  the  ideal 
N-wave  and  for  the  N-wave  as  reproduced  in  the  booth  using  predistortion. 

Using  these  functions  and  the  Johnson-Robinson  prediction  procedure,  the  com- 
pEirison  between  the  loudness  levels  of  both  signatures  was  made.  For  the  same 
rise  time,  overpressiire  and  duration,  the  difference  in  calculated  loudness 
was  found  to  be  about  0.4  phon.  Thus  the  residual  distortion  remaining  in 
the  reproduced  wave  - the  effect  of  iir5)erfect  cancellation  by  the  predistortion 
process  - is  seen  to  have  but  a small  effect  bn  siibjective  loudness. 

In  the  third  part  of  the  current  project  (see  Chapter  5-3),  hi^- 
pass  filtered  N-wave  signatures  were  \ased  in  order  to  explore  the  contribution 
of  the  low  frequencies  to  loudness . The  signals  were  generated  by  the  digital 
filtering  of  the  predistorted  N-wave  signature  in  the  computer.  In  Fig.  7 the 
full  process  of  generating  these  inpulses  is  shown  in  the  form  of  a simplified 
block  diagram.  The  N-wave  signature^  generated  in  digital  fbrm,  was  first 
transformed  into  the  frequency  domain  by  means  of  the  FFT  procedure  and  then 
predistorted  in  the  frequency  domain  (cf  above).  This  spectrum,  after  an 
inverse  FFT  into  the  time  domain,  was  used  as  the  reference  signed..  The  same 
predistorted  spectrum  (a  conplex  Fourier  transform)  was  alternatively  filtered 
by  a high  pass  digited.  filter  with  cut-off  frequency  of  either  25  or  50  Hz. 

This  spectrum,  after  an  inverse  FFT  into  the  time  domedn,  became  the  test 
signed. . 


The  digital  filter  simulated  a sliqple,  Ideed.  RC  circuit  (see  Fig.  12) 
with  the  following  frequency  response: 
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(3.9) 


H(.)  = ^ 


iJoilC 


where  RC  is  a time  constant  of  the  filter.  The  con^jlex  product  of  the  above 
function  and  the  Fourier  transform  of  the  predistorted  N-wave  signature  %((!>) 
gives  (after  converting  into  time  domain)  a digitally  calculated  signature 
representing  a signal  filtered  by  the  ideal  high  pass  filter, 


Fp(t)  =5='^{H(u)Fjj((o)) 


(3.10) 


The  resulting  pressure  signatures  recorded  inside  the  booth  are 
shown  in  Fig  13-  They  are  almost  identical  with  the  theoretically  calculated 
responses  of  the  high  pass  RC  filter  for  N-wave  signature  (see  Crocker  and 
Sutherland  [7^]) . Pulses  obtained  by  means  of  real  aneilog  filters  are  sli^tly 
different  because  of  the  non-ideal  characteristic  of  real  filters  (Hilton  and 
Newman  [75]).  In  Fig.  l4  the  comparison  between  the  calculated  and  analog 
pulses  is  shown. 

3 Measurements  of  Transient  Sounds 
3.^.1  Measurements  in  Time  Domain 


Measurements  of  the  pressure  time  history  of  the  transient  sounds 
imposes  special  conditions  on  the  equipment  used  for  this  purpose.  The  physical 
parameters  that  describe  lii^)ulslve  sounds  include:  pesik  sound  pressure  level, 
rise  time,  rate  of  decay,  pressure  variation  with  time  and  repetition  time. 

All  of  these  features  must  be  measured  and  recorded  accurately. 

Ideally  a transducer  and  a measuring  system  should  have  a linear 
frequency  response  between  DC  and  well  beyond  the  highest  significant  frequency 
in  the  signal.  However,  in  practice  measuring  equipment  always  has  limitations 
in  high  and  low  bands  of  the  frequency  response,  which  gives  as  a result  incorrect 
values  of  the  measured  parameters  of  the  Impulsive  signals. 

Tanlquchl  [76]  and  Hilton  et  al  [75]  have  defined  the  requirements 
for  the  i^er  and  lower  cut-off  frequency  of  the  measuring  system  as  the 
frequencies  where  the  Fourier  transform  f\aictlon  of  the  shock  pulse  drops  by 
40  dB. 


Crocker  and  Sutherland  [74]  defined  the  lower  and  vpper  limits  of 
the  equipment  for  measuring  real  sonic-booms  by  theoretical  calculations  of 
■the  filtering  effects  Vflpon  the  ideal  N-wave  signature  and  blast  waves,  without 
taking  into  account  the  absorption  of  high  frequencies  by  the  atmosphere  and 
other  aeabimtsina.  which  increase  the  rise  time  of  the  shock  wave.  They  recommended 
a linear  frequency  response  with  lower  and  upper  cut-off  frequencies  of  0.05  and 

10.000  Hz  respectively. 

In  the  current  e:q>erlment  the  range  of  the  frequency  spectrisn  for  all 
simulated  iiipiQ.ses  used  in  the  subjective  tests  extended  from  nearly  0 .1  to 

5.000  Hz;  thus  the  requirement  for  equipment  for  measuring  these  inpulses  was 
an  essentially  flat  frequency  response  covering  that  range. 
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The  neasui'ement  system  must  also  be  stable  and  sensitive  enough  to 
convert  the  pressure  variation  into  electrical  signals  which  can  be  detected 
or  stored  on  tape  without  degradation  or  contamination  by  noise.  The  trans- 
ducer and  electronic  equipment  should  have;  a good  phase  response,  less  than 
1.5  dB  ringing  and  should  be  completely  dan^jed  after  about  100  microsecond. 

The  measurements  should  be  made  in  a precisely  specified  environment  in  order 
that  the  experiments  coiold  be  duplicated  or  checked  by  other  researchers. 

Large  nearby  reflective  surfaces  should  be  avoided  and  the  angle  of  incidence 
of  the  transducers  should  be  kept  constant,  etc.  (Ross  et  al  [77]). 

All  measurements  in  oxu*  experiments  were  carried  out  at  a fixed 
position  of  the  microphone  (approximately  at  the  level  of  the  siibject's  ear) 
in  the  closed,  empty  booth.  The  microphone  enployed  v&b  a Bruel  and  Kjaer  1 
inch  condenser  microphone  type  4l46,  with  a random  incidence  corrector  type 
UA  00^^.  This  microphone  is  designed  for  sound  pressiu'e  measurements  down 
to  less  than  0.1  Hz  (this  is  achieved  by  closing  the  air  equalization  hole  of 
a standard  condenser  microphone) , and  the  random  incidence  corrector  makes  it 
practically  omnidirectional  ip  to  10  kHz.  Together  with  the  microphone,  the 
Bruel  and  Kjaer  carrier  system  type  263I  was  esployed.  The  frequency  range  of 
the  carrier  system  extends  from  DC  to  150  kHz,  and  it  provides  an  ideal  first 
stage  processing  of  sonic-boom  and  other  inpulses  that  carry  a strong,  very  low 
frequency  content. 

The  output  signals  from  the  carrier  were  either  recorded  back  on  the 
Briiel  and  Kjaer  FM  type  recorder  type  7003,  which  features  the  frequency  range 
at  speed  15  ips  from  DC  to  10  KHz,  or  photographed  by  a Polaroid  camera  from  a 
dual  beam  cathode-ray  Tektronix  osclllosccpe  with  storage  type  5103N.  This 
method  of  measuring  an  impulse  noise  involving  a microphone  and  oscilloscope 
is  regarded  as  very  acc\urate  and  is  widely  used  in  practice  (Martin  et  eJ.  [78]). 

3-4.2  Spectrum  Measurements  - Fast  Fom’ler  Transform 

In  order  to  calculate  the  loudness  of  transient  sounds  by  the  Johnson 
and  Robinson  procedure,  it  is  necessary  to  calculate  the  energy  spectrum  of  the 
waveforms.  This  was  done  in  the  computer  by  the  Fast  Foiurier  transform  procedure, 
which  was  also  used  for  generating  test  signals  (filtering  and  predistorting 
procedures) . 

A real  transient  inpulse  may  be  represented  as  y(t)  (real)  in  the 
time  domain,  or  equally  well  as  Y(u)  (complex)  in  the  frequency  domain.  Here 
Y((i))  is  the  Foturier  transform  given  by 

00 

Y(u)  = Jy(t)e"J"^  dt  (3.11) 

“00 

The  total  energy  of  the  inpulse  y(t)  is  finite  and  is  r^resented  by  the  time 
integral  of  an  instantaneous  "power" 


(3.12) 


From  Peurseval's  theorem. 
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y2(t)dt=i-  r (Y(u))  l^dw  (3.13) 

J ^ ” J 


-00  -09 

p 

The  quantity  |Y(u)J  is  called  the  energy  spectral  density;  and  when  plotted 
against  frequency  (positive  only)  gives  the  energy  spectrum.  The  area  under 
the  spectrum  curve,  multiplied  by  l/v,  must  equal  the  total  energy  W.  The 
dimensions  of  the  energy  spectral  density  are  energy  per  cycle  per  second. 

There  are  three  general  ways  to  evaluate  the  Fovrier  transform  of 
the  iii?)ulsive  sounds.  The  analytic  approach  involves  evaluating  the  integral 
(3 .11);  this  is  possible  only  when  the  analysed  waveform  is  represented  in 
the  form  of  an  analytical  equation.  In  case  of  the  experimentally  recorded 
signals,  it  is  possible  to  obtain  the  solution  only  for  simplified  analytical 
models  approximating  the  signature. 


The  second  approach  Involves  the  analog  analysis  of  recorded  signa- 
tures by  means  of  a set  of  analog  (real)  filters.  This  can  be  acccnplished 
by  recording  several  repetitions  of  the  analysed  pulse  on  a magnetic  tape  loop 
and  then  analysing  the  series  of  repetitions  by  a wave  analyser  (which  incor- 
porates the  filters;  Olsen  [79]). 

As  a third  approach,  the  analog  filter  techniqvie  is  replaced  by 
digital  analysis  in  computers  or  by  specially  designed  digital  spectral  analy- 
sers. using  the  procediire  that  we  have  referred  to  as  Fast  Fourier  transform 

(fft). 

In  general,  for  the  real-valued  record  y(t)  the  Fourier  transform 
integral  (3.13)  niay  be  restricted  to  the  finite  interval  (O,  T)  so 


X 

Y(f,  T)  = Jy(t) 


e-J^’^^dt 


(3.14) 


Assuming  that  y(t)  is  sanpled  at  "N"  equally  spaced  points  a distance 
"h”  apart  for  arbitrary  frequency  "f"  the  discrete  version  of  the  equation 
(3.i4)  can  be  shown  as 


Y(f,  T)  = h^y(nh)< 


-j27rfnh 


(3.15) 


t = nh,  n = 0,  1,  2,  . . .N-1,  T = Mh 


At  these  frequencies  the 
by 


k = 0,  1,  2,  N-1 

The  results  axe  unique  only  out  to  k = n/2,  the  Nyquist  cut-off  frequency. 

The  Fast  Fourier  transform  procediires  were  developed  to  con?)ute  the 
coefficients  of  the  equation  (3 .16)  with  reduced  ccn^uter  time  coo^ared  to 
standard  methods.  Most  often  used  procedures  are  based  on  the  Coolev-Tukev 
algorithm  {[80],  [8l],  [82]). 

Because  of  the  nature  of  the  digital  Fourier  transform  the  user  should 
be  aware  of  seme  limitations  and  errors  introduced  to  the  final  result:  these 
have  been  discussed  in  Chapter  3.3.1. 

In  the  practical  situation  only  a finite  block  of  data  can  be  processed 
by  the  FFT  procedure.  The  size  of  this  block  determines  the  frequency  resolution 
achieved,  of  step  size,  Af  = l/Nh  (where  Nh  = T,  total  length  of  the  san^le) . 

The  nuni>er  of  computed  Fourier  coefficients  'vrfien  operating  on  real  sampled  points 
is  exactly  one-half  the  nunfcer  of  sample  points  in  the  analysed  blodi  (u).  The 
^crease  of  the  sampling  rate  in  time  increases  the  rmyiTtpirn  frequency  range 
(Fmax  = l/2k);  but  for  the  same  frequency  resolution  the  computer  memory  and 
time  requirements  must  be  increased.  All  the  mentioned  factors,  i.e.,  cost, 
memory  size  available,  maximum  frequency  range,  sampling  rate,  frequency  reso- 
lution, etc.,  have  to  be  carefully  chosen  to  achieve  the  desired  result. 

The  sampling  rate  in  the  time  domain  is  usually  chosen  according  to 
the  required  frequency  range  of  the  spectrum.  The  basic  rule  is  that  the  sampling 
rate  must  be  at  least  twice  as  large  as  the  highest  frequency  component  in  the 
sampled  signal.  All  frequencies  beyond  that  limit,  if  existing  in  the  signal, 
should  be  filtered  out  because  of  the  danger  of  being  anslysed  as  lower  frequency 
components.  This  effect  is  called  aliasing. 

The  necessity  of  truncating  the  sampled  signal  in  the  time 
gives  rise  to  "leakage"  or  sidelobe  effects  in  the  frequency  domain,  i.e., 
spreading  of  the  main  lobe  of  the  true  power  spectral  density  function  cuid 
adding  an  Infinite  number  of  small  side  lobes.  Half  of  these  lobes  are  negative. 
Several  different  types  of  time  windows  are  used  Instead  of  the  simple  "box  car" 
function,  with  carefiiLly  designed  characteristics  for  smoothing  the  spectral 
samples  and  reducing  the  leakage  (see  Harris  [83]).  In  the  case  of  transient 
signals,  the  windowing  is  usually  made  with  the  simple  box  car  window  defined 
by 


transformed  V£d.ues  give  the  Fourier  componentB  defined 

Y(fj^,  T)  '^^Tin  I 

= {j =2^y(nh)  e (3.16) 

n=0 


u(t)  =0  t < -T 

* 1 -T  < t < T 
= 0 t > T 


24 


I 


All  spectrum  functions  of  the  test  pulses  in  our  e3q>eriments  were 
calculated  by  means  of  the  FPT  procedure  programmed  in  digital  computers 
(HP  2100  and  IBM  36O) . Some  typical  energy  spectral  density  functions  are 
shown  in  y^pendlx  6. 
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4.  EXPERIMENTAL  TECHMIftUE 

4.1  Design  of  Experiment 

4.1.1  General  Remarks 


The  objectives  of  any  psychological  experiment  are  to  measure  and 
interpret  psychological  events.  In  such  experiments  the  investigator  is 
lookuag  for  a relationship,  in  the  simplest  case,  between  two  variables,  the 
Independent  and  the  dependent  variable.  The  independent  variable  is  manipulated 
by  the  experimenter  and  the  dependent  one  is  manipulated  by  the  variations  of  the 
independent  variable. 

The  relationships  determined  by  the  experiment  must  be  also  appli- 
cable to  the  sitTxations  other  than  the  one  in  which  the  data  were  obtained; 
therefore  it  is  inportant  to  assure  the  generality  of  the  experiment.  The 
basis  for  the  generality  of  the  data  is  representative  sampling,  which  will 
provide  the  same  conposition  as  for  the  general  peculation.  There  are  several 
established  procedures  to  assure  representativeness  of  the  sample.  They 
include  random  sampling  with  replacement  method,  random  sampling  without 
replacement  (in  both,  each  siibject  from  the  popilLation  has  an  equal  chance 
of  being  selected) , and  rational  procedures  in  which  the  sample  is  composed 
in  the  same  way  as  the  general. population.  All  of  these  methods,  however, 
have  many  llmiitations  in  practice,  and  it  is  almost  never  possible  to  obtain 
truly  representative  samples  from  the  population. 

There  are  always  ntunerous  uncontrolled  variables  interfering  with 
the  investigated  relationships  in  the  psychological  experiments;  for  example, 
sex,  age,  elements  of  the  environment,  etc.  Keeping  the  experimiental  sample 
homogeneous  in  some  respects  (e.g.,  same  age  or  sex  in  the  experimental  groip) 
may  eliminate  the  extraneous  veiriables,  but  it  must  be  applied  with  caution, 
since  it  will  also  influence  the  generality  of  data. 

A good  way  of  minimizing  the  effect  of  uncontrolled  variables  is 
randomization.  This  is,  however,  most  effective  for  a large  number  of  samples. 
Nonselective  community  surveys  provide  an  example  of  this  approach. 

4.1.2  Loudness  Determination  Methods 

There  are  several  well  established  methods  for  loudness  determination 
available  in  experimiental  psychology.  Stevens  [84]  lists  seven  of  them.  The 
most  often  used  in  psychoacoustic  research  are; 

(i)  Method  of  adjustment  - observer  adjusts  stimulus  until  it 

is  subjectively  equal  to,  or  bears  some  desired  relation  to, 
the  criterion. 

(ii)  Minimal  change  procedure  - stimulus  is  varied  vpward  and/or 
downward  and  observer  communicates  its  relation  to  a criterion. 

(iii)  Paired  cottparlsan  or  constant  stimuli  procedure  - stimuli  are 
presented  in  pairs.  Each  stimulus  is  paired  with  another  or 
with  a reference.  Subjects  indicate  which  stimulus  in  each 
pair  is  greater  in  respect  to  a given  attribute. 
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(iv)  Rating  scale  or  magnitude  estimation  procedure  - each  of  a 
set  of  stimuli  is  given  an  "absolute"  nuaerical  rating  in 
terms  of  some  attribute. 

Each  of  these  techniques  has  specific  advantages  and  should  be  used 
according  to  the  purpose  of  the  investigation.  The  most  std table  for  loudness 
measurements  and  most  often  used  in  loudness  research  are  the  magnitude  estima- 
tion and  paired  con^jarison  methods.  They  are  frequently  combined  with  the 
adjustment  procedure. 

The  adj\istment  procedure  gives  the  subjects  control  over  the  stimuli 
throu^^i  some  sort  of  "transforming"  device.  Robinson  [85 J,  Stevens  and  Poulton 
[86],  and  Rood  (63],  among  others,  used  knobs  or  attenuators  as  a means  of 
control.  A controlling  system  regulated  by  a stbject  introduces  a bias  which 
will  be  dependent  ig)on  the  characteristics  of  the  device.  Stevens  and  Poulton, 
for  exanple,  reported  some  differences  between  results  obtained  with  a sone 
potentiometer  and  those  with  a decibel  attenuator.  The  use  of  the  adjustment 
procedure,  on  the  other  hand,  reduces  the  experimental  time  and  allows  more 
precise  measurement  of  loudness  levels. 

The  method  of  paired  comparisons  is  especially  suitable  as  a means  of 
establishing  subjective  equality.  It  gives  the  relative  levels  of  two  sounds 
at  which  they  are  considered  equally  loud  (or  equally  annoying)  . The  method 
is  considered  more  natural  than  the  others,  and  it  usually  gives  more  consistent 
results.  Yet  there  axe  some  specific  problems  associated  with  it.  Subjects 
tend  to  Judge  the  last  presented  stimulus  in  the  pair  as  less  acceptable  than 
the  first.  This  can  be  corrected,  however,  by  alternating  the  presentation 
order,  and  this  method  was  used,  for  example,  by  Shepherd  and  Sutherland  [6]. 

The  method  of  paired  conparisons  does  not  give  any  information  about 
quantitative  differences  between  two  compared  signals,  or  deviations  from 
equality  with  respect  to  some  criterion;  thus  it  is  usually  suitable  only  for 
relative  equality  level  measurements. 

The  magnitude  estimation  techniques  result  from  Stevens'  numerical 
scaling  procedures  [87],  [88].  Usually  the  standard  stimulus  is  presented  and 
observers  are  given  some  particular  number  describing  its  loudness,  noisiness 
or  annoyance,  or  they  may  be  free  to  choose  their  own  number.  Often  there  is 
no  standard  signal  in  the  experiment  and  the  observers  are  asked  to  assign 
nxaibers  which  appear  to  them  appropriate  in  proportion  to  a criterion.  It  has 
been  indicated  that  the  results  of  magnitude  estimation  tests  may  vary  with  the 
choice  of  an  intensity  of  a standard  stimulus  and  its  relation  to  the  numerical 
value  assigned  to  it  (Stevens  [88]  and  [86],  Heilman  and  Zwlslocki  [69]  and 

[90] ) , and  the  relative  frequency  with  which  the  standard  is  presented,  etc. 

Attempts  to  compare  the  techniques  of  paired  conparisons  and  magnitude 
estimation  have  indicated  that  both  methods  provide  highly  reliable  data.  While 
the  magnitude  estimation  technique  appears  to  be  the  faster,  the  paired  com- 
parison technique  seems  the  more  natural  for  the  subjects  (Clark  and  Kryter 

[91] ) • 

^.1.3  Statistical  Considerations 

An  eaqperlment  must  be  planned  so  that  the  resvilts  are  not  biased  by 
irrelevant  variables.  To  assvire  this,  the  investigator  should  Introduce 
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randosnization  into  the  experimental  procedvire  (e.g.,  in  order  of  presentation 
of  the  experimental  stimuli)  and  use  some  sorb  of  strictly  objective  measure 
of  error. 


There  are  well  known  methods  developed  in  parametric  statistics 
which  are  widely  used  for  testing  the  statistical  significance  of  psychological 
experiments.  They  allow  the  evaluation  of  \diether  the  independent  variable  has 
any  influence  on  the  dependent  variable  and  often  it  is  done  by  assuming  that 
there  is  no  such  influence.  This  notion  is  called  the  null  hypothesis. 

A freqixently  used  method  for  testing  the  null  hypothesis  is  a "t" 
test,  which  is  regarded  as  a most  powerful  test' in  statistics.  Two  general 
ass\mptlons  must  be  satisfied  for  usage  of  the  "t"  test.  The  sample  from 
which  the  data  aure  obtained  is  assvoned  to  be  normally  distributed  and  the 
variances  of  the  experimental  groups  eire  assumed  to  be  equal.  Recent  studies, 
however,  have  indicated  that  the  "t"  test  remains  accTirate  despite  violation 
of  these  assusptions,  providing  they  are  not  large  (Anderson  [92],  Boneau 
[93],  Sheridan  1^4]) . 

The  "t"  test  is  iisiially  applied  to  paired  sanples  (it  may  also  be 
applied  to  sanples  which  are  not  paired  - Govaden  [95]),  and  it  con?)ares  the 
standard  deviations  of  data  with  the  mean  difference  between  sauries  treated 
as  a VEirlate. 

Another  powerful  tool  in  statistics  used  routinely  in  experimental 
psychology  is  an  analysis  of  variance  and  Fisher  "F"  ratio  test.  The  analysis 
of  variance  is  also  based  on  the  assuaptlons  that  experimental  errors  are 
Independent  in  the  probability  sense,  have  equal  variance,  and  are  normally 
distributed.  The  conplete  analysis  of  variance  sorts  out  and  estiioates  the 
variance  components,  within  the  groups  and  between  the  groups,  and  in  the 
second  place  provides  for  a test  of  significance.  The  "F"  ratio  is  the  ratio 
between  ''between  groups"  variance  and  "within  groups"  Variance,  and  allows  one 
to  test  the  null  hypothesis  in  case  of  more  than  two  samples. 

4.2  Experimental  Procedure 

Our  investigative  work  consisted  of  four  experimental  test  series 
involving  human  observers.  Each  series  was  designed  to  test  the  subjective 
loudness  of  a different  type  of  transient  sound;  i.e.,  sonic-boom  N-wave 
signatures,  minimized  "low-boom"  signatures,  filtered  N-wave  slgnatwes  and 
idealized  "quarry  blast"  signatures. 

All  siibj active  tests  were  carried  on  in  the  UTIAS  Sonic  Boom 
Simulator  (for  more  details,  see  Chapter  3),  and  the  paired  conparison  tech- 
nique was  used.  All  test  signals  series  were  generated  by  conputer  and  then 
recorded  on  FM  tfpe  recorder;  thvis  the  subjects  taking  pcui;  in  the  comparison 
tests  had  no  control  over  the  compared  signals. 

During  the  13-20  minute  test  sessions  the  signals  In  pairs  were 
presented  to  each  suibject  while  seated  singly  in  the  simulator  booth.  They 
were  required  to  identify  which  sound  in  a pair  was  judged  to  be  the  louder, 
and  to  communicate  this  verbally,  through  the  intercom,  to  the  experimenter. 

Three  judgement  scores  were  used:  "louder",  "may  be  louder"  and 
"equal  loudness"  (similarly  for  annoyance  when  aiuioyance  weus  a criterion). 
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Thus  a set  of  five  nianerical  scores  was  obtained: 


Test  pulse  louder  = -2 
Test  pulse  may  be  louder  = -1 
Both  equally  loud  = 0 
Reference  pulse  may  be  louder  = 1 
Reference  pulse  louder  = 2 


The  reference  and  test  signals  in  pairs  were  presented  twice  to  the  stibjects, 
each  time  in  the  opposite  order;  i.e.,  AB  and  BA.  The  test  series  for  each 
coo^jared  pair  of  signatures  consisted  of  three  "warm  vq)"  pairs  at  the  beginning 
of  the  series,  and  several  cco5)arison  pairs,  each  with  different  overpressure 
ratios  between  the  test  and  reference  signals.  They  were  presented  to  the 
observers  in  random  order. 

The  test  series  were  generated  according  to  the  recommendations 
for  siibj active  comparison  tests  formulated  by  Reichardt  and  Niese  [9^].  The 
intervals  between  the  conpared  pulses  were  kept  between  500  and  1000  milli- 
seconds assuring  no  mutual  interference  and,  on  the  other  hand,  small  possible 
memory  errors.  The  duration  of  the  comparison  pulses  never  exceeded  2:1  ratio 
and  the  continuoizs  sound  in  the  last  part  of  the  experimental  project  was 
presented  for  400  ms. 

The  individual  results  for  each  value  of  the  tested  parameters  (e.g., 
rise  time  or  dilation)  were  plotted  in  the  form  of  graphs  - relative  loudness 
(in  scares)  versus  overpressure  ratio  between  the  test  and  the  reference 
signals.  Typical  exaaples  of  such  graphs  are  shown  in  Fig.  15.  From  each 
graph  the  oveipressxire  ratio  for  equal  loudness  (score  = 0)  was  determined 
and  the  average  of  these  values  over  all  popilations  of  participating  siibjects 
was  used  to  construct  the  final  equal  loudness  curves. 

ITEIAS  male  graduate  students  were  used  as  subjects  in  all  experiments. 
Audiograms  were  obtained  before  and  after  each  experimental  session  for  all 
observers.  In  addition  they  were  examined  prior  to  the  experiment  by  a qualified 
otolaryngologist  and  found  to  have  healthy  ears  and  hearing  levels  no  more  than 
15  dB  down  from  the  ISO  atidiometric  zero  [27],  except  for  occasional  dips  to 
-20  dB  in  the  frequency  range  above  kHz.  TMs  is  well  within  the  range  of 
hearing  (O  to  25  dB  dowi^  considered  normal.  Before  each  test  series  careful 
instructions  (see  Appendix  5)  were  given  to  the,t,individual  subjects. 

t 

All  the  enpirical  results  were  ccopared  with  the  theoretical  calcula- 
tions based  on  the  Johnson-Robinson  prediction  procedure  [2]  developed  for 
sonic-boom  signatures. 

The  analysis  of  variance,  "F"  Fisher  ratio  test,  and  "t"  test  were 
used  as  a standard  tool  for  testing  statistical  significance  of  the  experimental 
resixlts . 
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5-  SUBJECTIVE  EXFERIMEMTS 

5'1  Sitbjectlve  Loudness  of  N-Wave  Sonic  Boom 


5*1.1  Introduction 

There  is  continuing  interest  in  human  response  to  sonic-boom  type 
pressure  waveforms.  Ir.  particular , the  role  of  the  rise  time  and  duration  of 
the  N-wave  signatures  in  controlling  subjective  loudness  has  been  under  study. 

The  typical  values  of  scnic-boom  characteristic  parameters  measured 
at  ground  level  depend  upon  the  type  of  aircraft,  characteristics  of  the  flight 
and  atmospheric  conditions.  In  general,  the  peak  pressure  /p  varies  from  1 to 
10  psf  and  duration  varies  from  0.05  to  0.3  of  a second  (NASA  measurements  [75])* 
Typical  values  for  Concorde  SST  aircraft  are  1 to  2 psf  overpressure  and  from 
200  to  250  ms  duration.  Reported  rise  time  values  for  Concorde  range  from  0.1 
to  15  ms  [97].  Measurements  for  various  military  aircraft  made,  for  example, 
during  the  experiments  in  Edwards  AFB  appear  to  show  rise  time  in  a range  from 
0.2  to  45  ms.  True  rise  times  may  be  even  shorter,  as  0.1  to  0.2  ms  represent 
instrumental  rise  time  limits. 

Experimental  sxibjectlve  studies  with  sonic-boom  signatures  in 
laboratory  environments  have  been  conducted  by  Zepler  and  Harel  [1],  Lukas  and 
Kryter  [60],  Kryter  [65],  Shepherd  and  Sutherland  [6],  and  Rood  [63]. 

The  theoretical  prediction  procedure  for  determining  the  apparent 
loudness  of  sonic-booms  has  been  proposed  by  Johnson  and  Robinson  [2]  following 
earlier  proposals  by  von  Port  [39]>  and  Zepler  and  Harel  [1]  who  developed  their 
own  calculation  procedure  for  determining  the  loudness  of  sonic-boom. 


The  enpirical  studies  indicated  significant  change  in  loudness  of 
N-waves  with  variation  of  rise  time  and  marginal  change  in  loudness  with 
variation  of  duration.  Theoretical  predictions  based  on  the  spectrm  of 
sonic-boom  signatures  show  a similar  relationship. 

Shepherd  and  Sutherland,  as  well  as  Rood,  have  performed  tests  at 
different  reference  overpressure  levels.  Their  results  indicate  that  the 
level  appears  to  have  little  effect  on  the  relative  effects  of  rise  time  and 
duration  in  the  range  from  0.8  to  2.4  psf  in  Sh^herd’s  experiments  and  between 
0.5  and  2.0  psf  in  Rood's  tests. 


This  segment  of  the  present  investigation  is  very  similar  in  concept 
to  the  investigation  of  Shepherd  and  Sutherland.  However,  advances  in  the 
present  simulation  booth  design  allowed  for  a five-fold  shorter  rise  time. 


5.1.2  Results  and  Discussion 


! 


Two  separate  sequences  of  sonic-boom  conparlsons  featuring  N-wave 
signatures  were  carried  out  with  twenty  subjects.  Exanples  of  N-wave  booms 
r^roduced  in  the  Sonic  Boom  Simulation  Booth  and  used  in  the  cvirrent  experi- 
ment are  shown  in  Figs . 16  and  17) • 

In  the  first  sequence  of  experiments  the  boom  duration  "D"  was  held 
constant  at  200  ms,  the  rise  time  "t"  was  varied  over  the  range  0.22  to  10  ms, 
and  the  peak  overpressure  over  the  range  0.5  to  4 psf  (24  to  192  N/n^)  the  latter  only 
for  the  longer  rise  times.  For  each  rise  time  the  overpressure  of  the  test 
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N-wave  was  adjusted  so  that  the  observer  could  judge  for  which  overpressure 
level  the  loudness  matched  that  of  a reference  N-wave  with  1 ms  rise  time, 

1 psf  (48  N/h^)  overpressure,  and  200  ms  duration.  In  this  fashion  contours 
of  equal  loudness  versus  rise  time  were  developed. 

In  the  second  sequence  of  experiments  the  test  sonic-boom  rise  time 
was  held  constant  at  1 ms  and  its  dtu-ation  D was  varied  from  100  to  250  ms;  a 
second  equal  loudness  contoiir  (overpressure  ratio  vs  duration)  was  defined  by 
conqparison  tests,  adjusting  the  overpressures  from  0.5  to  2 psf  (24  to  96  N/m^) . 
The  reference  N-wave  was  the  same  as  the  previous  one. 

The  two  experimentally  determined  equal  loudness  contoxirs  for  the 
N-wave  signatures  are  plotted  in  Figs.  18  and  19 . The  overpressure  ratio  is 
defined  by  r^^io’ 

^level  ratio  ^°®10  p^^^  ^ ^°®10  1 psf 

Each  siibject  carried  out  approximately  I80  judgements  during  the 
course  of  two  test  sequences.  The  experimentally  determined  standard  deviation 
for  each  plotted  point  is  indicated  by  the  vertical  bars  on  the  graphs.  It  was 
noted  that  the  deviations  among  the  individual  conparison  results  increased 
progressively  as  the  differences  between  the  features  of  the  reference-boom 
signature  and  the  test-boom  signature  increased.  This  reflects  the  increased 
conparison  difficulties.  The  standard  deviation  is  typically  about  1 dB;  but 
for  booms  having  a diuration  of  250  ms  (rise  time  1 ms)  it  rises  tp  to  1.4  dB 
and  for  booms  having  rise  time  of  10  ms  (duration  200  ms)  it  is  about  3*3 

The  one-way  analysis  of  variances  performed  on  the  experimental  data 
along  with  the  individual  results  are  shown  in  Appendix  1.  Each  table  (i.e.. 
Tables  1 and  4)  represents  averaged  individual  results  over  the  AB  and  BA 
orders  of  presentation.  At  the  bottom  are  shown  the  calculated  mean  values 
for  all  subjects. 

The  analysis  of  variance  for  the  rise  time  tests  yields  a large  sum 
of  squares  representing  variability  among  the  population  means  (SSB),  but 
much  smaller  values  of  the  component  representing  variability  within  the 
samples  (SSE) , i.e.,  between  the  subjects. 

The  Fisher's  "F"  ratio  is  significant  at  the  significance  level 
0.01.  Therefore  the  null  hypothesis  can  be  rejected,  and  we  conclude  that 
the  experimental,  differences  in  the  rise  time  comparison  are  real  and  not  due 
to  chauice. 


In  the  case  of  the  duration  conparison  the  Fisher  "F"  ratio  for  data 
obtained  in  the  experiment,  excluding  the  sequence  with  250  ms  duration,  is 
much  smaller  than  the  value  required  for  rejecting  the  nxill  hypothesis.  There- 
fore the  null  hypothesis  cannot  be  rejected,  and  we  conclude  that  there  are  no 
real  differences  between  loudness  of  the  sonic  booms  over  the  duration  rauige 
100  to  200  ms.  The  same  Fisher  ratio  calculated  over  all  points,  including 
signatures  with  250  ms  dvuration.  Indicates  significant  difference  between  mean 
values  at  the  confidence  level  of  0.01. 
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Along  with  tbe  present  results  (labelled  Niediswieckl) , Fig.  18 
reproduces  the  experimental  resists  of  Shepherd  and  Sutherland  [6]  along  with 
the  calculated  theoretical  predictions  by  Zepler  and  Harel  [1]  and  by  Johnson 
and  Robinson  [2].  There  is  generally  good  agreement  over  the  common  range, 
essentially  within  the  error  bars.  The  predicted  decrease  in  loudness  with 
increasing  rise  time  is  very  marked  above  0.5  ms.  At  10  ms  rise  time  the 
results  are  scmewhat  divergent,  but  with  a large  experimental  \mcertalnty . 

In  the  rise  time  range  between  1 ms  and  10  ms  the  equal  loudness 
curve  falls  8 dB  which  conpetres  with  the  value  of  13  dB  in  the  Shepherd  and 
Sutherland  tests  and  9 dB  predicted  by  Johnson  and  Rbblnson's  theory.  In  the 
range  below  1 ms  the  equal  loudness  curve  rises  an  additional  3 dB  at  the 
lowest  rise  time  of  0.22  ms  investigated  in  the  present  test. 

In  the  range  of  durations  from  100  ms  to  200  ms  the  loudness  appears 
to  be  independent  of  duration  changes;  however,  for  durations  of  about  250  ms 
the  experimental  data  indicate  an  abrupt  rise  of  the  eqxial  loudness  curve. 

Both  Shepherd  and  Sutherland  (experiment)  and  Johnson  and  Robinson  (theory) 
find  a negligible  influence  of  the  dxaration  on  the  svbjectlve  loudness  for  the 
total  experimental  range. 

5.1.3  Conclusions 


The -present  results  appear  reasonably  consistent  with  earlier 
theoretical  and  experimental  subjective  boom  data,  except  for  the  effects  of 
the  longer  boom  durations  (in  excess  of  250  ms)  shown  in  Fig.  19.  The  substan- 
tial rise  in  the  equal  loudness  contour  in  this  case  remains  unexplained; 
however  we  may  speculate  that  this  effect  may  be  due  to  excitations  of  body 
vibration  by  the  boom.  Denpsey  and  Leathetwood  [98]  have  shown  experimentally 
that  the  human  body  sensitivity  for  -vibrations  (with  discomfort  as  a criterion) 
is  maxima]  in  the  frequency  range  from  3 to  8 Hz.  The  dominant  frequencies 
of  sonic  boom,  proportional  to  reciprocal  duration,  lie  in  this  range. 

The  essentially  good  agreement  with  previous  data  adds  confidence  to 
our  experimental  technique  as  well  as  to  the  existing  theoretical  methods  for 
predicting  the  subjective  loudness  of  N-wave  signatures  (especially  the 
Johnson  and  Robinson  procedure) . Compared  with  earlier  experience  these  are 
now  seen  to  validate  over  an  expanded  parameter  range  given  by 

Rise  time  0.22  to  10  ms 

Duration  100  to  250  ms 

5.2  Subjective  Loudness  of  "MLnindzed"  Sonic  Boom  Waveforms 
5.2.1  Introduction 

One  of  the  major  problems  that  has  limited  development  of  sxpersonic 
ci-vll  a-viation  is  annoyance  caused  by  sonic  booms  generated  by  overflights. 
Therefore,  a prominent  avenue  of  research  has  been  the  exploration  of  techniques 
for  minimizing  or  possibly  altering  the  diaracterlstlc  N-wave  produced  in  the 
f ar-field  by  present  day  sipersonic  aircraft . 

McLean  fTl  observed  that  for  sufficiently  long  aircraft  (e.g.,  > 300 
ft  (~  90  m))  the  "midfield"  sonic-boom  signature  may  not  have  evolved  fulTy 
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^ into  an  N-wave  at  ground  level.  He  suggested  that  the  shape  of  the  aircraft 

cross-section  and  lift  distribution  could  be  modified  to  optimize  this  midfield 
waveform  for  reduced  subjective  loudness. 

Hayes  [99]  pointed  out  that  in  the  real  atmosphere,  because  the  char- 
acteristics coalesce  more  slowly  than  in  the  uniform  atmosphere,  the  "midfield" 
signature  "freezes"  instead  of  evolving  into  an  N-wave  at  ground  level.  "If 
midfield  effects  persist  to  TiR/2  real- atmospheric  scale  heights  (H)  in  homo- 
geneous atmosphere,  they  will  persist  indefinitely  below  the  aircreift  in  real 
atmosphere."  Tailoring  of  the  midfield- type  sonic-boom  signature  for  minimum 
boom  was  developed  and  extended  by  George  and  Seebass  in  the  series  of  papers,  e.g., 

[8]>  for  flight  in  an  Isothermal  atmosphere.  They  ejqploited  the  phenomenon 
of  "freezing"  of  the  midfield  signatures. 

The  mathematical  theory  of  Seebass  and  George  has  been  extended  to 
apply  in  detail  to  the  real  atmosphere  by  Dcu:den  [9],  [10].  Her  formalism 
permits  minimization  of  either  the  initial  shock  of  the  signature  or  maximum 
overpressure,  compared  with  an  N-wave.  Again  this  is  acconplished  by  means  of 
an  especially  tailored  distribution  of  the  airersft  cross-section  and  lift. 

By  means  of  such  tailoring  Darden  computed  a fanily  of  minimized  or 
"low-boom"  signatures  associated  with  certain  proposed  " second  generation"  ' 

siQjersonic  transport  configuration  (Fig.  20)  The  expectation  was  that  for 
given  aircraft  volvnne,  wi^ight,  flight  altitude  and  Mach  number,  these  signatures 
should  soiaid  less  loud  than  norma]  N- waves.  Our  objective  in  the  present  test 
series  has  been  to  test  this  notion  experimentally.  The  signatures  would  be 
simulated  in  the  ITTLAS  Sonic  Boom  Simulator,  and  Jury  tests  of  the  siobjective 
loudness  would  be  conducted.  Finally  the  results  would  be  compared  with  cal- 
ciilated  values  by  the  Johnson  and  Robinson  method. 

Darden's  signatures  are  not  qviite  symmetric  fore  and  aft  (cf . Fig.  20)  . 

However,  the  Johnson-Robinson  method  for  predicting  loudness  is  predicated  on 
fore-aft  symmetry.  For  this  and  other  reasons  of  a practical  nature,  Darden's 
signatures  were  replaced  by  symmetric  ones  in  the  tests,  the  relationship  being 
as  in  Fig.  20.  The  differences  are  not  great  and  It  is  thoxight  their  effect  on  , 

the  subjective  loudness  of  these  "low-boom"  signatures  should  be  minimal. 

I- 

The  following  series  of  subjective  tests  attempt  to  establish  en^jiri-  i 

cally  the  relationships  between  the  subjective  loudness  and  "low-boom"  signature  j. 

characteristic  paxameters,  i.e.,  the  flat  top  duration  D^,  and  the  ratio  shock  , 

overpressure/peak  oveiT)ressure  (x  = Z!psh/^PMAX»  * i 

h 

r 5-2.2  Results  and  Discussion 

The  "low  boom"  or  "flat  top"  sonic  boom  signatures  did  not  reproduce 
accvirately  in  the  ITTIAS  simulation  booth,  despite  the  equalization  filter 
adjustments.  This  was  resolved  by  means  of  the  cco^uter-alded  scheme  for 
’ ' "predistorting"  the  input  signal  (see  Chapter  3-3  for  details).  The  resulting  , 

I ; typical  pressure  signatures  recorded  inside  the  simulation  booth  as  used  in  the 

subjective  tests  are  shown  In  Fig.  21. 

( In  all  test  sequences  the  signatures  bad  a fixed  rise  time  (1  ms)  and 

> ) total  duration  (150  ms) . 

r 
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Two  separate  cco^arlson  sessions  were  carried  out.  In  the  first  one 
the  flat  top  duration  of  the  slgnatxires  was  adso  held  constant  (Di  * 30  ms)  and 
the  ratio  x = ^sh/^PmaX  (^^o^t  shock  overpressure/maximum  overpressure  ratio) 
was  varied  within  the  range  0.2  to  1.0.  The  equal  loudness  contour  (overpressure 
ratio  of  test  and  reference  signals  versus  "x")  was  defined  through  the  coo^jarl- 
sons  of  these  signatures  with  the  reference  N-wave  signature  having  the  same  rise 
time  (t  = 1 ms)  and  total  duration  (D  = 150  ms),  and  overpressure  of  = 0.5 
psf  (24  N/n^) . Ten  observers,  all  UELAS  male  graduate  students,  took  part  In 
this  experiment;  each  of  them  carried  out  about  120  judgements. 

In  the  second  test  sequence  the  overpressure  ratio  x = ^sh/^MAX 
was  held  constant  at  the  level  x = 0.5  and  the  equal-loudness  contour  of 
overpressiire  ratio  vs  flat  top  duration  (D^)  was  determined  for  the  "low-boom" 
signatures  having  a flat  top  duration  within  the  range  10  to  60  ms  at  a total 
duration  of  150  ms  (l.e.,  from  0.0667  to  0.4  of  the  total  duration).  The  refer- 
ence N-wave  had  the  same  total  duration  (150  ms)  and  rise  time  (l  ms)  as  previously 
but  the  overpressure  was  fixed  at  1 psf  (48  N/n^) . Eight  observers  took  part  In 
this  experiment,  and  each  carried  out  about  100  judgements  during  the  test  sequence 

Two  equal-loudness  contours  derived  from  the  experimental  results  for 
the  "low-boom"  signatures  are  lll\istrated  In  Figs.  22  and  23.  The  first  shows 
the  overpressure  level  ratio  vs  x = i^Psa/^MAX*  second  shows  the  over- 

pressure level  ratio  vs  flat  top  duration  D^.  The  overpressure  level  ratio  Is 
defined  by 


^level  ratio  ^°®10 

where  = overpressure  of  the  reference  N-wave, 

= maximum  overpressure  of  the  test  "low-boom"  signature. 

The  plotted  eqxial-loudness  curves  are  based  on  the  averaged  values 
calculated  from  the  e:q)erlmental  results  for  each  subject.  The  vertical  bars 
Indicate  the  experimentally  determined  standard  deviation.  The  standard  devia- 
tions for  the  ratio  "x"  conparlsons  are  within  the  range  0.8  to  1.5  dB;  for  the 
flat  top  dviration  conparlson  they  range  from  0.6  to  1.4  dB.  The  detadled 
results  for  each  experimental  series  and  all  individual  results  are  shown  in 
Appendix  2,  along  with  the  result  of  the  one-way  analysis  of  variance. 

The  Fisher  ratio  "F"  is  significant  for  variable  "x"  cooparisons 
at  the  significance  level  0.01;  thus  the  differences  between  the  means  for  the 
different  values  of  ratio  "x"  are  taken  to  be  real. 

The  analysis  of  variance  performed  on  the  results  of  the  second  test 
series  (with  varying  flat  top  dwatlon  Di)  does  not  permit  rejection  of  the  null 
hypothesis;  thus  the  differences  between  mean  values  of  the  subjective  loudness 
for  different  flat  tcp  durations  are  taken  to  be  the  result  of  chance. 

It  was  fo\uid  that  to  malntedn  equal  loudness  the  overpressure  ratio 
must  be  Increased  by  11.7  dB  as  the  parameter  "x"  Increases  from  0.2  to  1.0 
(Fig.  22) . The  actxial  properly  sceiled  waveforms  judged  as  equally  loud  are 
shown  In  Fig.  24.  The  conparlson  suggests  that  subjective  loudness  of  the 
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"low-boom"  type  of  signature  depends  mainly  i^jon  the  front  (and  rear)  shocks. 
This  conclusion  is  considered  valid  as  long  as  the  ratio  t/Di  is  of  the  same 
order  as  used  in  the  e:q)eriment  (l.e.,  t/Di  =0.033) » or  smaller. 

The  results  (for  the  particular  value  of  ratio  rise  time  and 

duration  used  in  the  test)  are  well  approximated  by  the  formula: 

(N-wave)  = ZlPgjj  + 0.11  ("low-boom")  (5-3) 

for  equal  loudness,  provided  the  N-wave  has  the  same  rise  time  and  duration 
(Fig.  25).  This  tells  us  that  the  peak  pressure  iiPMAX  contributes  only  one 
ninth  as  much  to  the  loudness  as  the  front  shock  (and  similarly  for  the  rear 
shock);  that  is,  the  front  (and  rear)  shock  aii5>litudes  (for  fixed  rise  time) 
do  dominate  the  loudness.  The  influence  of  the  subjective  loud- 

ness will  be  higher  for  higher  values  of  ratio 

For  the  effect  of  the  "flat  top"  duration  on  the  subjective  loudness 
refer  back  to  Fig.  23.  The  overpressure  level  ratio  for  equal  loudness  varies 
less  them  1 dB  with  increase  in  the  flat  top  duration  fromi).0667  to  0.4  of  the 
total  duration.  This  change  is  within  the  range  of  the  experimental  error  and 
was  not  found  statistically  significant.  Therefore,  we  can  infer  that  the 
duration  of  the  flat  top  part  of  the  "low-boom"  signature  has  a negligible 
Influence  on  the  subjective  loudness  in  the  experimental  range  of  values 
(similarly  for  the  ratio  t/Di) . 

The  results  of  the  "low-boom"  cooparison  tests  were  supported  by  the 
theoretical  loudness  calculations.  The  loudness  of  each  signature  judged 
equally  as  loud  as  the  reference  N-wave  was  calculated  from  the  energy  spectrim 
density  function  obtained  through  the  FFT  procedure  in  a digital  computer.  A 
few  exarples  of  the  energy  spectrum  for  the  "low-boom"  signatures  used  in  the 
experiment  are  shown  in  Appendix  7. 

The  Johnson-Robinson  method  for  N-wave  sonic-booms,  in  a version 
based  on^he  Stevens  Mark  VI  procedure  for  continuous  sounds,  was  followed  in 
the  calcinations . The  loudness  was  calculated  for  the  positive  parts  of  the 
signatures  only,  which  is  justified  after  Johnson  and  Robinson  on  the  ground 
that  the  'separation  between  front  and  rear  shocks  is  sufficiently  long  conpared 
to  the  crltic8Ll  time  of  the  auditory  syst^. 

The  results  of  these  calculations  are  conpared  with  the  calculated 
loudness  of  the  reference  N-waves  in  Figs.  26  and  27.  The  calculated  loudness 
(in  phons)  for  all  studied  "low-bocaa"  signatures  differs  from  the  calculated 
loudness  of  the  reference  N-waves  which  sound  equaQ.ly  loud  by  less  than  1 phon. 
This  very  good  agreement  of  the  enpirlcal  and  theoretical  results,  in  terms  of 
the  relative  lotidness,  supports  the  viability  of  the  Johnson-Robinson  loudness 
conparisons  between  N-wave  and  the  "low-boom"  family  of  signatures  within  the 
range  of  parameters  given  by: 

0.0667  < (D./D  = Dp/D)  <0.4 

~ at  D = 150  ms,  T = 1 ms 

0.2  < (x  * ^gjj/^j4Ax^  £ 1.0 
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A series  of  jury  tests  of  the  perceived  loudness  of  the  "low-boom" 


sonic-boom  signatures  has  been  conpleted  and  the  resvtlts  coinpaured  with 
theoretical  predictions.  The  caii5)arisons  indicate  that  the  loudness  of  these 
signatures  is  dominated  by  the  aji5)litude  ^SH  front  (and  rear)  shocks. 

The  peak  ajiplitude  can  thus  be  much  larger  than  that  of  an  N-wave  that  sounds 
equally  loud  (assuming  that  the  ratio  t/Di  is  of  the  same  order  as  in  the 
experiment  or  less) . Put  another  way,  an  N-wave  of  the  same  peak  aaplitude 
will  sound  much  louder  than  the  "low-boom"  signatures  with  small  ratio  "x". 

Based  on  Darden's  calculations  [9],  [10]  of  possible  "low-boom"  signatvtres 
for  realizable  aircraft,  with  ^sh/^^PMAX  half,  the  attainable 

loijdness  reductions  are  rovighly  equivalent  to  those  restalting  from  halving 
the  present  N-wave  anplitudes. 

There  is  an  inportant  caveat  concerning  the  above  results:  they 
refer  solely  to  sonic  booms  as  heard  outdoors.  Indoor  sonic  booms  are  quite 
different;  .their  waveforms  as  well  as  anplltudes  have  been  grossly  modified  by 
the  transmission  characteristics  of  walls  and  windows.  These  are  normally 
dominated  by  low  frequency  resonances,  and  are  quite  insensitive  to  the  shock 
aaplitudes  and  rise  tiaies.  The  major  parameters  governing  Indoor  boom  intensity 
are  the  iiipid.se  (area  under  positive  half  of  outdoor  waveform)  and  duration.  For 
the  "low  boom"  signatures  these  parameters  are  not  much  altered  conpared  with 
the  standard  N-wave;  indeed,  the  iapulse  may  be  Increased.  Thus  the  "low  boom" 
signatures,  although  much  quieter  outdoors  than  N-waves  for  ccoparable  size 
aircraft,  offer  no  advantage  Indoors. 

The  relative  loudness  predictions  of  the  Johnson-Robinson  theory 
conformed  very  well  to  the  measurements.  Thus  their  potential  for  applica- 
bility to  a much  broader  range  of  transient  sounds  is  indicated. 

The  above  results  make  it  clear  that  the  rise  time  of  the  shock  waves 
in  inpulsive  sounds  is  a major  parameter  controlling  the  subjective  loudness, 
along  with  their  peak  anplitude. 


Subjective  Loudness  of  Filtered  N-Wave  Signatxire 
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5. 3*1  Introduction 

Analysis  of  the  energy  spectral  density  of  a typical  N-wave  slgnatiue 
reveals  that  there  are  two  specific  break  points  in  the  gross  envelope  of  the 
spectrum  function  (see  Fig.  1,  after  Johnson  euid  Robinson  [2]) . The  envelope, 
starting  from  the  low  frequency  end  of  the  spectrum,  rises  at  6 dB  per  octave 
vp  to  a peak  freqxiency  of  Hz  (D  = duration),  then  feJLls  at  6 dB  per 

octave  and  finally  from  a frequency  of  l/irr  Hz  (t  = rise  time)  drops  at  12  dB 
per  octave.  Typically,  the  first  maximum  occurs  in  the  region  ranging  from 
2 Hz  to  6 Hz  for  current  svpersonlc  aircraft  (i.e.,  dtiratlon  ranges  from  100 
to  300  ms) . Thus  most  of  the  energy  of  the  sonic-boom  N-wave  signature  is 
concentrated  far  below  the  normal  audio  frequency  range  of  the  human  ear.  We 
know  from  well  established  evidence  that  the  sensitivity  of  the  ear  is  extremely 
low  at  frequencies  below  the  range  of  20  Hz,  called  "Infrasonlc" . Whittle  et  al 
[28]  have  shown  that  the  decrease  of  the  equal-loudness  curves  in  the  range  from 
2^  to  3 Hz  is  approximately  60  dB  at  the  levels  from  threshold  \p  to  70  phons. 

Experimental  laboratory  investigations  of  human  and  animal  response 
to  sonic-boom  are  usvtally  carried  out  in  ^eclally  designed  loudspeaker  facili- 
ties in  order  to  fujJLy  reproduce  N-wave  signatirres,  including  the  Infrasonlc 
content.  Thus  frequency-wise,  it  is  necessary  to  extend  the  flat  response  of 
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! such  a simulation  facility  down  to  neaxly  DC.  This  can  be  done  by  defeating 

the  usual  low-frequency  roll-off  by  use  of  either  an  airtight  chamber  or 
earphones  with  an  effective  pressxire  seal  (see  Chapter  3).  The  qt;iestion  arises; 
What  influence  upon  the  subjective  human  response  has  the  intense  infrasonic 
part  of  the  sonic-boom  spectrxmi  in  the  situation  of  whole  body  immersion  in 
the  acoustic  field?  If  this  Influence  were  found  to  be  small  or  negligible,  it 
would  be  possible  to  avoid  bvdlding  the  costly  simulating  equipment;  instead 
the  subjective  sonic-boom  tests  could  be  conducted  with  relatively  inexpensive 
high  fidelity  aiodio  eqvdpment  with  frequency  response  from  20  to  10,000  Hz. 

An  indication  that  the  energy  below  approximately  40  Hz  has  an 
insignificant  effect  vqjon  the  loudness  was  given  by  Zepler  and  Harel  [1]  in 
1965.  They  used  an  analog  (real)  high  pass  filter  with  cut-off  frequencies 
ranging  from  20  to  240  Hz  and  reproduced  the  filtered  N-wave  signatiares  with 
specially  designed  earphones.  Introduction  of  high  pass  filters  with  cut-off 
frequency  of  either  20  or  40  Hz  was  found  to  have  no  effect  on  the  subjective 
lovidness  of  sonic  booms  of  rise  times  of  3 and  1 ms.  The  results  of  Zepler  and 
Harel  agree  with  later  theoretical  estimations  by  Rood  [63].  He  calculated 
the  correlation  between  the  energy  in  discrete  bands  of  the  sonic  boom  spectrum 
and  the  subjective  loudness  level  of  the  stimvili.  Correlations  in  the  frequency 
bands  of  5-10  Hz,  10-20  Hz  and  20-40  Hz  were  foiand  to  be  insignificant;  the  40-80 
Hz  band  was  fovind  significant  at  the  5%  level,  and  the  remainder  of  the  level 
below  1%.  The  range  of  parameters  used  in  this  compaxison  were;  rise  time 
varying  from  2 to  16  ms  and  peak  pressure  varying  from  0.5  to  2 psf. 

The  present  test  series  is  essentially  a repetition  and  refinement 
of  the  Zepler-Harel  experiment  [l]  in  the  situation  of  whole-body  immersion  in 
the  acoustic  field.  The  analog  high  pass  filters  axe  replaced  by  digital 
filtering,  and  the  filtered  waveform  is  further  ingprcrved  by  the  "predistortion" 
technique.  The  resiolts  on  subjective  loudness  axe  compared  with  the  predictions 
of  the  Johnson-Robinson  theory  [2]  to  provide  a further  test  of  that  theory. 

5*3*2  Results  and  Discussion 

Fifteen  male  subjects  were  tested  with  loudness  as  a Judgement  criterion 
and  ten  with  annoyance  as  a Judgement  criterion.  All  of  them  had  been  checked, 
as  described  earlier,  for  normal  hearing,  and  all  were  UTIAS  graduate  students. 
Each  subject  carried  out  about  4o  JuuJgements. 

The  paired  comparison  technique,  as  in  the  previous  experiments,  was 
used.  The  reference  N-wave  with  1 psf  (48  H/n^)  overpressure,  1 ms  rise  time 
and  150  ms  diiration  was  generated  in  "predistorted"  form  by  a computer.  Two 
test  signatures  were  obtained  from  the  pre&istorted  reference  signature  by 
means  of  a digital  high  pass  filter  having  two  cut-off  frequiencies;  25  and  50 
Hz.  The  details  of  the  filtering  process  have  been  given  in  Chapter  3*3. 

Typical  pressure  signatuires  recorded  inside  the  booth  axe  shown  in  Fig.  l4. 


In  the  first  test  series  the  overpressuire  ratio  levels  defined  as 


^level  ratio  ~ ^°®10  ^ 


= 20  log, 


(5*4) 


* ) 


^ref  ^ 

were  found  for  the  conditions  of  equal  loudness  for  the  respective  filtered 
signatures  compared  with  the  reference  N-wave.  Based  on  the  values  of  the 
respective  signature  overpressures  Juidged  as  equally  loud,  the  absolute 
subjective  lovidness  levels  for  these  signatvtres  were  calcvdated  using  the 
Johnson-Robinson  [2]  procedure.  Both  experimental  and  theoretical  results  are 
shown  in  Table  3*1*  The  first  row  contains  the  empirically  established  over- 
pressures of  the  signatvires  Judged  by  the  observers  as  equally  loud  (values 
shown  represent  the  averaged  results) . Comparison  indicates  that  the  loss  in 
infrasonic  frequency  content  due  to  filtering  must  be  compensated  by  an  Increase 
in  overpressure  for  equal  loudness.  The  differences,  which  increase  with  cut-off 
frequency,  are  however  very  small.  ^7 


standard  deviations  calcvilated  from  the  test  data  were  0.88  dB  for 
a cut-off  frequency  of  25  Hz  and  I.36  dB  for  a cut-off  frequency  of  50  Hz. 

This  reflects  growing  difficulties  in  loudness  con5)arison  for  higher  cut-off 
frequencies  because  of  the  increasing  differences  in  character  between  con?)ared 
sounds. 


The  indlvidiaal  results  along  with  statistical  analysis  are  shown  in 
Appendix  3*  The  "t"  test  applied  to  the  results  of  the  loudness  cooparisons 
show  no  significant  difference  between  mean  values  obtained,  indicating  no  real 
difference  in  loudness  between  filtered  N-wave  signatures  for  the  two  cut-off 
frequencies  of  the  hl^  pass  filter:  25  and  50  Hz. 

The  theoretica31y  calculated  loudness  levels  of  the  reference  so\md 
and  test  signatures  (second  row  in  Table  3*1) > which  were  found  to  be  equally 
loud  through  the  subjective  test,  differ  by  less  than  0.6  phon.  This  good 
agreement  provides  further  svqiport  for  the  viability  of  the  Johnson-Robinson 
loudness  prediction  procedure,  especially  their  special  method  for  the  evalua- 
tion of  the  low-frequency  bands  (see  Chapter  2.4.3). 

All  svibjects  dviring  the  loudness  conparison  session  were  asked  about 
the  difference  in  character  between  the  test  and  reference  signals  in  terms 
other  than  loudness  parameters.  All  remarked  on  relatively  easily  detected 
differences  in  "sharpness",  frequency  content,  and  subjective  duration.  This 
led  us  to  the  conclusion  that,  even  though  the  subjective  loudness  of  sonic 
boom  does  not  depend  upon  the  energy  below  50  Hz  range,  other  subjective  para- 
meters may  be  affected  by  this  part  of  the  spectrum. 

Another  test  sequence  similar  to  the  one  described  above  was  carried 
for  ten  siibjects  with  annoyance  as  a criterion  of  judgement.  The  overpressures 
for  equal  annoyance,  averaged  over  all  observers  (third  row  in  Table  3*1) j show 
an  opposite  tendency  from  the  case  of  loudness  conparison.  To  state  it  simply, 
annoyance  was  slightly  increased  as  a result  of  the  filtering  out  of  the  low 
frequency  sound  energy. 

A similar  phenomenon  has  been  reported  by  Turner  and  Bums  [100] 
who  foiaid  that  svbjects  comparing  a 5 second  long  burst  of  one-third  octave 
band  white  noise  centred  at  5150  Hz  and  the  same  burst  with  additional  one- 
third  octave  band  white  noise  centred  at  the  much  lower  frequency  of  250  or 
500  Hz  jvidged  the  first  one  as  the  more  annoying. 

The  standard  deviation  for  the  annoyance  comparison  test  was  for 
both  cases  equal  to  1.0  dB. 

The  individual  results  are  shown  also  in  appendix  3*  In  the  same 
appendix  are  shown  the  results  of  the  statistical  analysis  performed  on  the 
eaperlmental  data.  The  "t"  test  indicates  no  significant  differences  between 
means  at  cut-off  frequency  of  25  and  50  Hz  for  annoyance  as  a judgement 
criterion.  The  same  "t"  test  performed  over  the  results  obtained  in  the 
tests  with  different  criteria  (i.e.,  loudness  and  annoyance)  but  the  same 
cut-off  frequency  does  not  allw  rejection  of  the  null  hypothesis  for  the 
cut-off  frequency  of  25  Hz,  but  gives  the  t ratio  valiie  close  to  the  value 
required  for  the  0.1  level  of  significance.  For  the  cut-off  frequency  of 
50  Hz  the  calculated  value  of  "t"  allows  rejection  of  the  null  hypothesis  at 
the  0,1  level  of  significance  and  is  very  close  in  value  to  the  "t"  parameter 
required  at  the  0.05  level  of  significance. 
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Therefore,  the  differences  between  the  resvilts  of  the  subjective 
e^qperiments  with  loudness  and  annoyance  as  Judgement  criterion  are  acc^ted 
as  real  and  dvie  to  factors  other  than  chance  at  a 0.1  level  of  significance, 
at  a cut-off  frequency  of  50  Hz. 

5.3.3  Conclixsions 

A series  of  test  coii?)axisons  between  unflltered  and  high-pass  filtered 
N-wave  signatures  (cut-off  frequencies  of  25  and  50  Hz)  was  carried  out.  The 
loudness  level  of  the  filtered  signatures,  with  their  reduced  low  frequency 
content,  was  found  to  be  only  marginally  lower  than  that  of  the  comparable 
vmflltered  N-wave  (both  signatures  had  a rise  time  of  1 ms,  duration  of  150  ms) . 
This  result  was  consistent  with  the  theoretical  loudness  predictions  based  on 
the  Johnson-Roblnson  [2]  method.  It  was  also  found  cospatlble  with  the  previous 
resvilts  obtained  by  Zepler  and  Harel  [1]  (experimental)  and  Rood  [63]  (calculated). 

The  above  result  conformed  very  well  to  one  of  the  general  conclvislons 
of  the  N-wave  investigation  in  Chapter  5*1:  there  was  no  significant  change  in 
subjective  loudness  with  variation  of  the  N-wave  duration  (except  for  the  anomalous 
abr\;q>t  rise  in  excess  of  250  ms) . Increased  duration  increases  the  energy  in  the 
infrasonic  frequency  range:  thus  insensitivity  to  duration  translates  into  an 
insensitivity  to  the  infrasonic  frequencies  in  the  N-wave. 

The  second  test  sequence  with  annoyance  as  a criterion  indicated  a 
slightly  increasing  tendency  in  annoyance  with  reduction  of  the  low  frequency 
energy  in  the  N-wave  signatiure.  A siMlar  phenomenon  was  r^orted  by  Turner 
and  Bums  [100] . 

The  difference  between  loudness  and  annoyance  of  the  filtered  signa- 
tvires  were  found  statistically  significant  at  0.1  level  of  significance  for  the 
case  of  high  pass  filtering  with  cut-off  frequency  of  50  Hz. 

5 Subjective  Loudness  of  Idealized  "Blast"  Signatures 

5.^.1  Introduction 

In  three  preceding  paragraphs  (i-e.,  5.1»  5*2  and  5*3)  it  has  been 
shown  that  the  relative  loudness  predictions  based  on  the  Johnson  and  Robinson 
[2]  theory,  agree  very  well  with  subjective  loudness  measurements  for  three 
types  of  transient  sounds;  the  N-wave  signature,  the  "low  boom"  signature  and 
filtered  N-wave  signatirre. 

The  last  part  of  this  project  was  designed  to  extend  the  applica- 
bility of  this  theory  to  even  broader  types  of  transient  sounds  and  also  to 
verify  the  absolvcte  values  of  the  theoretically  calculated  loudness  levels  by 
the  Johnson  and  Robinson  method. 

Considerable  attention  has  beer,  devoted  recently  to  research  on  quarry 
blasts'  effects  on  human  population,  by  various  government  agencies,  due  to  the 
nisnerous  conplalnts.  For  this  reason  the  Impulsive  sounds  chosen  for  subjective 
testing  were  based  on  the  quarry  blasts,  measured  for  example  by  Taylor  et  al 
[11].  The  range  of  the  blast  signature  parameters  recorded  in  Taylor's  project 
extends  for  peak  overpressure  from  0.42  psf  (20.2  N/n?)  to  4 psf  (192  N/m^)  and 
for  duration  of  the  signature  from  0.l4  to  2.7  seconds.  An  example  of  pressiure- 
timie  history  given  in  this  report  is  shown  in  Fig.  28. 
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The  pressiire  signatures  of  a series  of  quarry  blasts  show  great 
individual  variations.  Flg;ire  28  is,  in  fact,  more  regular  in  its  features 
than  most.  Thus  it  has  been  thoxight  worthwhile  to  suppress  the  variation  in 
signature  shape  by  idealizing  the  waveform.  The  idealized  signature,  shown 
in  Fig.  29,  resenibles  the  real  signat\are  of  Fig.  28  in  having  a series  of 
shocks  or  echoes  (so  that  it  looks  like  a comb);  however,  the  "teeth"  of  the 
comb  are  all  of  identical  shape  and  spacing.  More  precisely,  the  electrical 
input  signatures,  generated  in  the  computer  without  "predistortion",  consist 
of  a different  nuniber  (M)  of  successive  pulses  each  25  ms  long  with  shock- 
like rise  time  of  0.22  ms  and  quasi -exponent! ally  decaying  "tail",  according 
to  the  formula: 


V = (1  - t/D)  e"®*^ 


(5.5) 


where  D = duration  of  the  Individual  pulse  (25  ms), 
a = constant  (O.l  ms~^) . 

The  input  signatures  together  with  their  spectral  density  fiuictions 
are  shown  in  Appendix  7,  and  the  pressure  waveforms  recorded  within  the  booth, 
and  used  for  the  subjective  tests,  are  shown  in  Fig.  29. 

5.^.2  Results  and  Discussion 

In  the  series  of  siibjective  cco^arison  tests  (using  techniques 
described  in  Chapters  3 and  4)  the  rise  time  of  each  individual  pulse  was  held 
constant  at  0.22  ms  and  the  duration  of  the  pulse  sequence  was  varied  from 
25  to  400  ms;  this  was  done  in  five  st^s,  the  five  different  signatures  being 
sequences  of  1,  2,  4,  8 and  16  indivldiml  pulses,  re^ectively.  The  peak 
overpressure  was  varied  over  the  range  0.35  psf  (I6.8  N/m^)  to  1.5  psf  (72  N/h^) . 

The  reference  pulse  in  each  comparison  series  was  always  another  blast 
signature,  being  twice  as  long  as  the  test  one.  Thus  four  con^jEirison  tests  were 
carried  out  and  the  duration  (t)  of  the  signals  in  pairs  were: 

25  ms  vs  50  ms 

50  ms  vs  100  ms 

100  ms  vs  200  ms 

200  ms  vs  400  ms 

This  arrangement  eOLlowed  us  to  keep  the  duration  ratio  between  con^ared  signals 
within  the  recommended  value  2:1  (see  Relchard  and  Niese  [96])»  since  a much 
larger  disparity  in  duration  makes  loudness  coo^arlsons  difficvilt.  Furthermore, 
this  procedure  permitted  covering  a wide  range  of  durations. 

In  this  fashion  contours  of  equal  loudness  versus  dxuration  (or  number 
of  individual  pulses  in  the  signature)  were  developed.  Fifteen  subjects  took 
part  in  the  series  of  tests,  each  making  about  90  cooparlsons.  The  individual 
results  along  with  their  statistical  analysis  cure  shown  in  Appendix  4.  The 
results  as  shown  again  represent  the  overpressure  ratio  of  reference  and  test 
signature  for  equal  loxidness,  averaged  over  the  AB  and  BA  order  of  presentation. 


The  resvilts  averaged  over  the  total  popvilation  of  the  observers 
were  then  transformed  Into  the  ratios  between  each  tested  pulse  and  the  new 
reference  level  - being  the  overpressure  level  of  the  longest  ("continuous") 
signature,  with  duration  of  400  ms.  The  overpressure  ratio  is  defined  by 
^level  ratio*  where 


^level  ratio  " -2°  losio  M » 1,  2,  4,  8 (5.6) 

The  resxilting  experimentally  determined  equal  loixlness  curve  (overpressure 
level  ratio  vs  duration)  is  plotted  in  Fig.  30  re  400  ms  case  taken  as  110 
phons  ( solid  line) . 

The  standard  deviation  varies  from  0.75  dB  for  ccn^iarisons  between 
pulse  sequences  of  200  ms  and  400  ms  duration,  to  1.45  dB  for  con5>arisons 
between  pulse  sequences  of  100  ms  and  200  ms  duration. 

The  Fisher  "F"  ratio  and  analysis  of  variance  performed  on  the  data 
(shown  in  Appendix  4)  allows  us  to  reject  the  null  hypothesis  at  0.01  level  of 
significance  and  to  conclude  that  the  differences  between  mean  values  obtained 
from  the  duration  comparison  are  real. 

The  overpressure  level  ratio  versus  duration  contour  rises  about  2 dB 
for  each  doubling  of  duration  in  the  range  25  to  100  ms.  For  signatures  longer 
than  100  ms  the  increase  is  progressively  slower,  approaching  as  an  asymptote 
the  level  for  continuous  sound  (in  excess  of  400  ms) . This  type  of  relationship 
between  the  loudness  and  the  duration  of  the  inpulsive  sounds  has  been  reported 
previously  in  the  literature  and  reflects  the  tenq^oral  integration  effect  of  the 
auditory  system  (Munson  r^5]j  Garner  and  Miller  [40],  Green  et  al  [4l],  Harris 
[42],  Ploop  and  Bouman  [43],  etc).  It  was  discussed  in  Chapter  2.3.2. 

As  mentioned  before  (see  Chapter  2.3.2),  there  is  little  agreement 
in  the  literature  as  to  the  value  of  the  time  constant  of  the  human  auditory 
system.  The  estimates  range  from  20  to  200  ms  and  there  are  conflicting 
reports  whether  the  temporal  integration  function  depends  \q)on  the  spectrum  of 
the  stimuli  and  its  intensity  level.  It  is  rather  hard  to  accurately  measure 
the  critical  time  of  the  ear  for  the  averaged  data  obtained  in  the  current 
test,  because  of  the  relatively  large  intervals  between  the  experimental 
values  of  duration;  however  we  may  estimate  it  as  being  approximately  between 
100  and  200  ms. 

The  observation  that  the  critical  time  for  the  test  "blast"  pulses, 
estimated  from  the  empirical  equal-loudness  curve,  is  longer  than  100  ms  has 
a (pacific  consequence  for  the  theoretical  prediction  of  the  absolute  loudness 
levels  based  on  the  Johnson-Robinson  method.  They  assumed  [2]  that  the  time 
necessary  to  evoke  full  response  of  the  human  auditory  system  is  70  ms  - the 
value  estimated  by  von  Port  [39]*  The  underestimation  of  the  real  vdue  of 
the  crltlceLL  time,  as  used  in  calculations,  should  lead  to  a significant  over- 
estimation  of  the  absolute  values  of  loudness  levels  ced-culated  by  the  Johnson- 
Robinson  scheme.  Indeed,  the  calculated  loudness  levels  for  the  test  "blast" 
signatures  Judged  as  equally  loud,  using  time  constant  value  equal  to  70  ms 
are  overestimated  by  approximately  4 phons  ccnpared  to  the  experimental  cvurve 
(Fig.  30  ■ curve  3)»  in  terms  of  the  absolute  loudness  levels.  However,  in 


terms  of  the  relative  values,  the  predicted  levels  are  in  very  good  agreement 
vd.th  the  data  acquired  in  the  experiment  up  to  tc  = 70  ms:  this  can  he  seen 
by  downshifting  the  70  ms  curve  by  3-6  phons,  which  superimposes  it  on  the 
experimental  curve.  The  absolute  loudness  scale  in  phons  on  the  graph  is 
set  up  by  the  calculated  loudness  level  of  the  signature  (using  Stevens'  Mark 
VI  procedure)  with  the  duration  of  400  ms,  which  can  be  regarded  as  a "contin- 
uous sound"  (see,  for  example,  Reichard  and  Niese  [96])  for  human  auditory 
system  response.  (The  term  "effective  loudness"  is  used  on  the  graphs  as  a 
reminder  that  the  experimental  curve  is  still  in  terms  cf  the  Zp  level  ratio 
for  equal  loudness,  as  defined  by  equation  5*6.) 

A similar  result  (i.e.,  overestimation  of  the  absolute  siibjective 
loudness  by  the  original  Johnson -Robinson  theory)  was  reported  by  Rood  [63], 
who  compared  the  experimenteJ.  measurements  of  the  absolute  values  of  siibjective 
loudness  of  half  N-wave  signatures  with  predictions  by  the  Johnson-Robinson 
theory.  The  graph  showing  this  comparison  is  reproduced  in  Fig.  31.  The 
theoretical  calculations  show  constant  overestimation  of  the  loudness  by  about 
3.2  dB,  which  is  very  close  to  the  result  obtained  in  this  investigation. 

Agreement  in  both  absolute  and  relative  loudness  levels  between 
the  experimiental  data  and  the  Johnson-Robinson  predictions  is  attained  when  we 
change  the  critical  timie  value  in  the  calculations  to  175  nis  (this  lies  in  the 
range  estimated  previously  from  our  experimental  data);  cf.  Fig.  30  again, 
curve  2. 


The  Johnson-Robinson  principle  of  linear  integration  of  sound  energy 
with  timie  applies  only  for  signal  duration  T ^ t^.  For  T > tc  it  is  implicit 
that  the  judged  loudness  is  constant  at  the  value  for  continuous  sound  of  the 
same  spectral  content.  This  accounts  for  the  loudness  "ceiling"  of  110  phons 
shown  in  Fig.  30  for  the  Johnson-Robinson  calculation  (curve  2) . The  linear 
integration  assumption,  it  is  seen,  staxts  to  break  down  when  the  duration 
exceeds  two  thirds  of  the  critical  time  Tq. 

By  contrast  the  experimental  loudness  curve  (solid  line  in  Fig.  30) 
increases  smoothly  with  duration  toward  the  continuous  sound  asymptote.  It 
is  as  though  there  were  some  sort  of  saturation  effect  in  the  neurological 
response.  Without  attempting  any  specific  modelling  inferred  from  observed 
cochlear  and  neurological  mechanisms,  one  can  postulate  a simple  e:9onential 
approach  to  satxiration.  Thus,  in  place  of  Johnson  and  Robinson's  [2]  linear 
law  which  defines: 

0)2 

band  pressure  level  = ^J'  L dw  x T > 

we  introduce  the  generalized  definition 


The  factor  in  brackets  is,  in  effect,  an  interpolation  function:  it  approaches 
T/Tq  at  small  signed,  durations,  and  imlty  at  large  signal  d\irations.  There  is 
precedent  for  forms  similar  to  (5*8)  in  earlier  work  (e.g.,  Munson  [45],  Plonp 
and  Bouman  [43],  Hughes  [101],  etc). 
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• The  Johnson-Robinson  method  [2],  as  generalized  to  incorporate  equa- 
tion (5*8) > has  been  applied  to  the  simulated  blast  waves  deedt  with  in  Fig. 
30.  Several  choices  of  integration  time  were  tried,  and  beat  fit  to 
experiment  was  found  with  Tq  = 125  ms.  This  is  labeled  ciorve  3 in  Fig.  32. 

For  caii5)arison,  curve  2 is  a calculation  by  the  original  Johnson-Robinson 
method,  best  fit  being  obtained  by  the  choice  Tq  = 175  ms. 

In  Fig.  33  the  solid  curve  ( eiqjerimental)  and  curve  3 (generalized 
Johnson-Robinson  theory)  are  repeated  from  Fig.  32.  Curve  2 on  this  figure  is 
a calculation  by  the  original  Johnson-Robinson  method  with  the  choice  Tq  = 125 
ms.  This  curve  2 is  clearly  not  a best  fit  to  e:q>erlment,  on  an  absolute 
basis.  But  the  curve  is  a good  fit  on  a relative  basis:  if  it  is  downshifted 
by  1.3  dB  it  matches  the  solid  curve  well  below  T = 100  ms.  In  fact  the  down- 
shifted Tq  = 125  ms  curve  2 of  Fig.  33  is  virtually  identical  with  the  best- 
fit  Tg  = 175  ms  curve  2 of  Fig.  32. 

We  have  seen  that  the  assun5)tlons  = 125  ms  and  Tq  = 175  ms  both 
give  conqjarably  good  relative  fits  to  esqieriment  (below  their  respective  upper 
limits) ; this  underlines  the difficulty  of  inferring  the  correct  physical 
critical  integration  time  by  such  a matching  procedure  based  on  the  linear 
integration  assimgption.  The  exponential  integration  assunptlon  (5*8),  with 
Tc  reinterpreted  as  an  e“l  time  constant,  appears  to  yield  a much  more  precise 
determination . 

5.^.3  Conclusions 


A set  of  equal-loudness  cooparison  tests  for  idealized  "quarry 
blast"  waves  has  been  carried  out.  The  waveforms  were  simulated  as  a sequence 
of  shock-decay  pul8esJ\|\|\J\j  the  duration  (~nunt)er  of  shocks)  being  varied. 
The  cxurve  of  overpressure  ratio  in  phons  vs  duration  shows  a steep  rise  at 
small  duration  (<  100  ms,  inpiU.sive  sound),  but  the  curve  bends  over  to  approach 
an  asynptotic  vpper  limit  at  long  duration  (>  400  ms,  continuous  sound). 


The  Johnson-Robinson  predictive^ theory  for  inpulsive  sound  [2]  gives 
good  absolute  agreement  with  the  low  end  of  the  duive^ -provided  the  critical 
integration  time  of  the  ear  tq  is  changed  from  70  ms  (their  -i^-ue)^  ,tp  175  ms. 
Values  of  tq  smaller  than  175  ms  give  good  relative  agreement,  but  the^vtrve- 
is  shifted  vpward.  In  this  theory  the  band  pressure  levels  are  assumed  to 
increase  linearly  with  signal  duration  T via  a factor  (T/tq) . 


In  ovu*  proposed  generalization  of  the  Johnson-Robinson  theory  the 
band  pressvure  levels  are  assumed  to  follow  an  exponential  saturation  law 
(1  - e"™/’'’o) . The  generalized  theory,  with  this  redefined  tq  chosen  as  125  ms, 
shows  a very  good  fit  to  the  enplrlcal  curve  over  the  whole  range.  Thus  it 
appears  - at  least  for  the  "quarry  blast"  waveform  - to  Integrate  loudness 
prediction  methods  for  transient  and  continuous  sounds. 

The  question  arises,  can  this  generalized  Johnson-Robinson  theory  be 
applied  in  hindsight  to  the  earlier  tests  herein;  will  it  correctly  predict  the 
loudness  of  sonic-boom  wave-forms  of  the  three  kinds  tested?  The  answer  is 
yes,  since  the  generalized  and  original  Johnson-Robinson  theories  (with  different 
choices  for  tc  (see  abovej)  agree  fairly  well  at  the  smaller  valixes  of  T/tq 
appropriate  to  sonic  boom. 
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6 . CONCLUDING  REMAEUCS 


The  findings  and  conclusions  drawn  with  respect  to  each  of  the 
four  test  series  have  been  discussed  in  detail  at  the  ends  of  the  respective 
chapters.  A summarized  account  is  given  verbally  in  the  opening  Abstract. 

An  additional  overview  is  given  in  schematic  form  in  Table  6.1;  the  four 
waveforms  tested  are  sketched  at  the  head  of  four  coliimnsj  and  the  respective 
key  findings  and  conclusions  are  summarized  underneath. 

The  resid-ts  in  each  series  were  coopared  with  theoretical  predictions 
by  the  method  of  Johnson  and  Robinson.  All  but  the  long-diuration  quarry  blast 
judgements  were  foiond  to  be  in  very  good  agreement  in  terms  of  relative  loudness 
levels.  With  an  ad  hoc  "exponential  saturation"  modification  (including  adjust- 
ment of  the  critical  integration  time  of  the  ear)  the  predictive  method  was 
extended  to  enconpass  the  long  duration  signals  as  well.  Thus  the  applicability 
of  the  method  has  been  demonstrated  for  other  types  of  transient  soimds  than 
the  N-wave;  and  the  extension  to  the  method  appears  to  bridge  the  range  between 
inpulslve  and  continuous  soxands  of  similar  spectral  content. 

As  stated  earlier,  the  res\alts  for  all  loudness  judgements  in  this 
project  were  acquired  with  an  all  male  graduate  student  population  (which 
inposes  a limitation  on  age  of  about  21  to  30  years)  . In  order  to  generalize 
the  results  we  must  know  what  is  the  influence  of  age  and  sex  upon  the  auditory 
response  to  the  transient  soiands. 

The  reports  in  the  literature  indicate  that  the  effect  of  aging  upon 
the  auditory  system  (i.e.,  presbycusis  and  sociocusis)  is  negligible  below 
^ kHz  and  is  not  of  major  significance  below  5 kHz  (Robinson  and  Dadson  [104], 
Pollack  [105] , etc) . The  spectrum  of  test  pulses  in  our  investigation  was 
regarded  as  significant  only  up  to  5 kHz;  consequently  the  generalization  of 
our  data  over  a wider  age  group  shoiild  not  introduce  a serious  error. 

Some  investigators  have  found  small  differences  between  female  and 
male  listeners  in  aural  sensitivity,  especially  in  the  high  frequency  range 
(e.g..  Rood  [63],  Bauer  [I06]) , while  others  like  Robinson  and  Dadson  [104] 
have  fovuid  no  systematic  effects  of  sex.  Thus  the  error  introduced  by  using 
only  male  observers  in  our  limited  frequency  range  should  be  small;  however, 
more  research  in  this  direction  is  necessary. 
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Table  6.1  Subjective  CoBparlaon  Testa  - Sumnary  of  Restilts 
[Definition:  Oveirreesure  lerol  ratio  = -20  logio(  ^test/'^reference^  equal  loudness  ^ 
H-Wave  i Flat-Top  Boom  j Filtered  N-Wave  I Idealized 
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APPENDIX  1 

EXPERIMEOTAL  RESULTS  AND  STATISTICAL  ANALYSIS 


N-WAVE  COMPARISON  TESTS 
Table  1 

Results  of  Sonic-Boom  Comparison  Experiment.  Rise  Time  Variation 


Reference  SiKnature:  Duration  - 200  ms;  Rise  Time  - 1 ms 


Rise  time 
(ms) 


Level  Difference  (dB)  Required  for  Judged  Equality 


Observer  No. 


k 1 4 2 3.3  3 

3 1 2 2.7  2 2.5 

0 0.2  0 -1  0.5  0 

-4.6  -3.8  -5  -4.5  -6  -5-3 

10  -6  -5  -5.5  -11  -6 


Rise  time 
(ms) 


Rise  time 
(ms) 


Overall  Average 
All  Observers 
(dB) 


2.9 

2.4 

-0.09 

-4.7 

-7.8 


Level  Difference  (+dB)  Reqxiired  for  Judged  Equality 


Observer  No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0.2 

-11 

-11.25 

-10.75 

-11 

-9.1 

-11 

-11.75 

-11.25 

-9.25 

-9.5 

0.4 

-7 

-5.25 

-6.5 

-5.25 

-4.75 

-6.75 

-5 

-5.5 

-4.75 

-6.5 

0.6 

-3.5 

-1.25 

-4 

-5 

0 

-4 

-2 

-4.6 

-3 

-3 

0.8 

-1 

1.35 

-0.4 

0 

-1.25 

-2 

-0.75 

-0.4 

-0.8 

-1 

1.0 

-0.5 

3 

1 

0 

1.5 

0 

2.5 

0 

3 

2 

X 


Overall  Average 
All  Observers 


(dB) 


0.2 

0.4 


0.6 


0.8 


1.0 


-10.48 

-5.71 

-3.03 

-0.46 

1.26 


Observer  No. 

Level  Difference  (+dB)  Reqioired  for  Judged  Equality 

Ref.  D = 50  ms 
Test  D = 25  ms 

Ref.  D = 100  ms 
Test  D = 50  ms 

Ref.  D = 200  ms 
Test  D = 100  ms 

Ref.  D = 400  ms 
Test  D = 200  ms 

1 

2 

3 

3 

1 

2 

1.3 

4 

3 

1 

3 

1.35 

2 

1 

0.75 

4 

1 

1 

0 

0.33 

5 

2 

2.65 

0 

-1 

6 

2 

2 

1.25 

0.75 

7 

4 

3 

-1 

0.65 

8 

1 

1 

4 

2.75 

0 

9 

1 

3 

3 

2 

10 

1 

1 

0 

0.65 

n 

1.2 

2.35 

1.35 

2 

12 

3 

2 

2 

1.25 

13 

1 

3.35 

2 

1 

l4 

0 

2.75 

-1 

0 

15 

2 

2 

2.8 

0.33 

Overall 

Average 

1.59 

2.54 

1.34 

0.71 

Table  2 


Standard  Deviation  for  Data 
from  Blast  Signature  Conmarison  Experiment 


ms)  - D^  (ms) 

Standard 

Deviation 

(dB) 

50  - 25 

0.96 

100  - 50 

0.91 

200  - 100 

1.46 

400  - 200 

0.76 

Table  3 

One-Way  Analysis  of  Variance  for  Blast  Signattire  Duration  Comparison 


Re  400  ms  case  taken  as  0 dB 


Source  of 
Variation 

Degree  of 
Freedom 

Sxan  of 

Square 

F 

Between 

samples 

3 

SSB=267.77 

22.14 

Error 

56 

SSE=225.75 

Total 

59 

SST=493.52 

-2 


APPENDIX  ^ 


INSTRUCTION  AND  CONSENT  FORM  FOR  EXPERIMENIATION 
UTIAS  SONIC  BOOM  PROJECT 


"CONSENT" 


^consent  to  take 

part  in  a research  project  being  carried  out  at  the  Itaiversity  of  Toronto,  Institute 
for  Aerospace  Studies,  under  the  direction  of  Professor  H.  S.  Ribner.  I imderstand 
that  my  participation  will  inclxide  the  following; 

(a)  I will  be  exposed  to  pairs  of  transient  pressure  signals  (simiilated  sonic  booms 
and  other  similar  pulses)  while  seated  in  the  UTIAS  sonic  boom  simulation  booth. 

(b)  I will  ccmpare  the  relative  loudness  (or  annoyance)  of  the  paired  transient 
pressure  signals. 

(c)  Each  session  will  be  no  longer  than  15  minutes  and  there  will  be  no  more  than 
two  sessions  per  day,  totalling  up  to  six  sessions  per  test. 

( d)  A qxialified  otolaryngologist  will  examine  my  hearing  before  the  experiment  and 
in  addition  an  audiometric  test  will  be  performed  before  and  after  each  session. 

I understand  that  I may  withdraw  from  the  project  at  any  time. 


Dated  at  Toronto  this  day  of 


197  . 


Witness 


Signature 


INSTRUCTIONS 


The  p\irpose  of  the  test  is  primarily  to  determine  the  loudness  of 
different  transient  sovmds  as  a function  of  seme  particular  features  of  those 
pulses. 

You  will  hear  a series  of  ia^ulsive  sounds  while  seated  inside  the 
booth.  The  sounds  will  occiir  in  "pairs"  and  will  be  separated  by  few  seconds. 
Yo\ir  task  will  be  to  Judge  which  sound  in  each  pair  is  louder.  After  each 
pair  of  signals  there  will  be  an  interval  of  silence  so  you  will  have  enough 
time  to  communicate  your  Judgement  through  intercom. 

Please  base  your  decisions  on  your  subjective  feeling  of  loudness, 
disregarding  other  differences  between  pulses  which  may  occur.  In  order  to 
help  you  in  your  task,  five  standard  scores  will  be  used: 

first  (in  pair)  louder 
first  (in  pair)  may  be  louder 
both  equal 

second  (in  pair)  louder 
second  (in  pair)  may  be  louder 

Please  concentrate  on  the  Judgement.  There  is  no  "good"  or  "bad" 
answer;  we  are  Just  interested  in  your  siibjective  feeling.  Take  this  experiment 
seriously,  otherwise  it  will  be  a waste  of  time  for  you  as  well  as  for  us. 
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APPENDIX  6 

EXAMPLES  OF  COMPUTER- CALCULATED  SPECTRAL  ENERGY  DENSITIES 


OF  TEST  SIGNALS 
Part  2 

Spectral  energy  density  (calculated  by  a coii5)uter)  of  the  ideal 
filtered  N-wave  signat\ires  used  in  the  subjective  test. 
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APPENDIX  6 

EXAMPLES  OF  COMPUTER-CALCULATED  SPECTRAL  ENERGY  DENSITIES 


OF  TEST  SIGNALS 


Spectral  energy  density  (calciilated  by  a coniputer)  of  the  recorded 
booth  responses  to  ideal  "blast"  signatures  used  in  the  subjective 
test. 
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APPENDIX  7 
COMPUTER  FROGRAie 
Paxt  1 

Listing  of  the  FFT  Fortran  pregram  (IBM  360)  calculating  decibel  values 
of  one-third  octave  band  pressure  levels  using  the  Johnson-Robinson 
method  [2].  Program  reads  input  data  (signatures  stored  in  digital  form 
by  the  HP  2100  cco^juter  with  A/D  converter)  from  a magnetic  tape,  calculates 
Fourier  transform  using  a standard  FFT  subprogrEuu  and  confutes  energy  in 
one-third  octave  bands.  These  energies  are  converted  into  pressure  levels 
(using  the  Johnson-Robinson  scheme)  and  printed  as  output  data. 
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APPENDIX  8 


JOHNSON  AND  ROBINSON  METHOD  FOR  CALCULATING  LOUDNESS  OF  IMPULSIVE  SOUNDS; 

STEP-BY-STEP  COMPUTER  PROCEDURE 


1.  The  measured  iii?)ulsive  so\iRds  recorded  on  an  FM  tape  recorder  are  converted 
to  digital  form  by  an  A/D  converter  and  stored  in  a magnetic  tape  storage 
(HP  2100A) . The  frequency  of  samplinjg  is  fg  = 20  kHz.  The  maximum  amplitude 
is  scaled  to  ± 204?  discrete  points. 


2.  The  Fast  Fourier  svibroutine  CQn5)utes  the  cooplex  discrete  values  of  Fourier 
components 'for  k = 0,  1,  N-1  defined  as: 


where 


N-1 


exp 


n=0 


-J 


2iilsn 

N 


(1) 


x^  = x(n  . At) 


X(fjj,T) 

~~3 


T = total  duration  of  sample 
At  = sampling  time 

N = n\md)er  of  discrete  points  in  sample 
^ = ££,  frequency  s^aration  in  frequency  domain 

3.  The  energy  in  l/3  octave  bands  is  calculated  by  sinmlng  up  all  discrete 
values  of  the  energy  spectral  density  |Xit  • Ab  multiplied  by  the  frequency 
separation  A*)  = Ptt  • 

B 

\ ‘il  (2) 

k“A 

where  A and  B are  discrete  limits  of  a l/3  octave  band. 
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4 . The  effective  l/3  octave  band  pressure  levels  are  confuted  from  band 
energies: 


10  log 


10  2 
^ref 


(3) 


\diere 


= 0.0002  pbar 


t 


c 


I 


original  Jdhnson-Robinson  procedure 


generalized  Johnson-Rcbinson  procedure 


- effective  duration  of  the  signeil 

C 

T - crltlced.  time 
c 

In  the  case  of  sonic  booms  with  duration  > 100  ms  the  front  and  rear  shocks 
are  noimally  heard  separately  ("bang-bang") , “For  this  clrcvonstance  Johnson  and 
Robinson  reduce  the  presstire  levels  in  l/3  octave  bands  by  3 dB: 


= Lp  - 3 as  (I.) 

This  correction  is  restricted  to  sonic  boams>  or  double  In^ulses  of  similar 
duration. 


6.  Levels  in  bands  below  50  Hz  are  reduced  to  levels  at  this  frequency  (i.e., 
50  Hz)  which  produce  the  same  effect  (same  loudness  in  phons);  equal  lotidness 
ciurves  are  xised,  e.g.,  ISO  Recommendation  R226  [27].  These  weighted  levels  are 
then  combined  with  that  already  existing  in  the  50  Hz  band,  using  normal  rules 
of  decibel  addition. 


7 . The  resulting  modified  effective  pressure  levels  in  l/3  octave  bands  are 
lused  for  evaluating  the  svibjective  loudness.  The  standard  loudness  calculation 
procedure  by  Stevens  (ISO  Recommdnation  R532  [30])  was  used  in  the  present  work. 
The  loudness  in  sons  (St)  is  converted  to  loudness  level  in  phons  (P)  by  the 
equation 

P = 40  + 10  logg  S^  (5) 
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SONIC  BOOM  SIMULATION  FACILITY 

chamber,  (B)  Electronic  network  and  eqviipment 
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riG . 4 SIMPLIFIED  BLOCK  DIAGRAM  OF  EQUALIZING  NETWORK  FOR  LOUDSPEAKER- 
DRIVEN  SONIC  BOOM  SIMULATOR  BOOTH. 


DIGITAL  COMPUTER 


I 


SIMULATOR 


ntDICAIES  ADDITIONAL  DIGITAL  FILT  BRING  SUBROUTINE  USED  IN  THIRD  PART  OF  PROJECT 


f 


FIG.  8 EFFECT  OF  "rREDISTOBIION"  OF  ELECTRICAL  INPUT  SIGNALS  (TOP)  TO  UTIAS 
SONIC  BOOM  IN  ACHIEVING  DESIRED  WAVEFORMS  (BOTTOM)  RECORDED  BY  MICRO- 
PHONE IN  BOOTH. 

(a)  "Low-boon"  sonic  boom:  t = 1 ms,  D = 150  ms,  Di  = 30  iiis,  x = 0.8 
(B)  N-wave:  t = 1 ms,  D = 150  ms 


FIG.  9 EXAMPLES  OF  PRESSURE  SIGNALS  RECORDED  BY  MICROPHONE  IN  UTIAS 

SIMULATION  BOOTH  WITHOUT  (TOP)  AND  WITH  (BOTTOM)  PREDISTORTION 
OF  INPUT  SIGNAL. 

(a)  N-wave:  t = 1 ms;  D = 150  ms;  /'ip  = 1 psf 

(B)  "Low-boom"  sonic-boom:  t = 1 ms;  D = 150  ms;  D^  = 30  ms; 

X =1.0;  APmflv  = 1 psf. 
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FIG . 12  SIMPLE  RC  HIGH  PASS  FILTER  SIMULfflED  IN  DIGITAL  FORM 
IN  COMPUTER. 


A. 

fc=O.I  Hz 


fc  = 25  Hz 


C. 

f c » 50  Hz 


FIG.  13  EFFECT  CF  HIGH  PASS  FILTERING  ON  SONIC  BOOM  SIGNAL  (PRESSURE 
WAVEFCRIC  RECORDED  IN  UTIAS  SIMUL/fflION  BOOTH).  N-WAVE 
SIGNATURE:  t 1 ms;  D > 1^0  ms;  ^ « 1 psf. 


(a)  No  filtering  - effective  fg  * 0.1  Hz 

(B)  Filtered  - fc  = 25  Hz 

(C)  Filtered  - f » ■ 50  Hz 


A) 


B) 


FIG.  l4  COMPARISON  BETWEEN  EXPERIMENTAL  (ANALOG)  (A)  AND  FREDHTTED  (b)  TIME 

HISTORIES  FOR  SONIC  BOOM  MEASURED  BY  A SYSTEM  WITH  FREQUENCY  RESPONSE 
FROM  10  Hz  to  10  kHz  (FROM  REF.  74)  . 


FIG . Ip  RELATIVE  LOUDNESS  SCORES  VS  OVERPRESSURE  RATIO  BETWEEN  REFERENCE  AND 
TEST  BOOM. 


Siibject  No.  4 - duration 
Sid>Ject  No.  4 - duration 


200  ms;  rise  time  » 1 ms 
200  ms;  rise  time  - 0.^  ms 


AP  ItMl  ratio 
[dB] 


MEOZWIECKI  exp. 

SHEPHERD.  SUTHERLAND  expi 
JOHNSON,  NOBMSON  theo. 
ZEPLER,  HAREL  theo 


riM  tiRM  IimI 

FIG.  l8  EQUAL  LOUDMESS  CURVE.  TRADE-OFF  BETWEEN  OVERPRESSURE  lEVEL  RATIO 
(-20  logio  ^teet/^ref)  AND  RISE  TIME  FOR  200  ms  DURATION  N-WAVES 

^ ^ N/n?),  T = 1 ms. 


AP  level  ratio 


FIG.  19  EQUAL  LOUDNESS  CURVE.  TRADE-OFF  BEIEWEEN  OVERPRESSURE  LEVEL 
RATIO  (-20  logio  ^te8t/4>ref)  AND  DURATION  FOR  1 ms  RISE 
TIME  N-WAVE. 

/p  = 1 psf  (48  N/m^),  D = 200  ms. 


RT  - RiM  Tim* 

D|  - Front  Flat  Top  Duration 


PRESSURE  D - Duration 


FIG.  20  Tmi!AT.T7.m  vs  AITAINABLE  "LOW-BOOM”  SONIC-BOOM  SIGNATURES. 

"ATTAINABLE"  SIGNIFIES  REALIZABLE  VIA  AIRCRAFT  DESIGN  AND 
FLIGHT  PROCEDURE. 


FIG.  21  REPRODUCTION  CP  "LOW-BOOM"  SIGNATURES  BY  UTIAS  SIMULATION 
BOOTH  AS  USED  IN  TESTS  (WUH  PREDISTOKriON)  . RISE  TIME 
B 1 ms,  DURATION  = 1^0  ms. 

(a)  Ratio  front  shock  overpressvire/maxlmum  overpressure 
>0.^;  flat  top  duration  ■ 10  ms 

(B)  Ratio  front  shock  overpressure/maximum  overpressure 
B 0.^;  flat  top  duration  « 60  ms 

(C)  Ratio  front  shock  overpressure/maxlmum  overpressure 
> 1.0;  flat  top  duration  > 30  ns 
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(phons) 


A 


N’WAVE 


— AP  = 0.5  psf 
RT  s I ms 
D s 150  ms 


FIG.  26  CALCULATED  LOUDNESS  CF  SIGNATURES  JUDGED  EQUALLY  LOUD.  SOLID  LINE; 

LOUDNESS  OF  LOW-BOOM  SIGNATURES  VS  RATIO  FRONT  SHOCK  OVERPRESSURE/ 
MAXIMUM  OVERPRESSURE. 


LOUDNESS 


DI=D2 


FIG.  27  CALCULATED  LOUDNESS  OF  SIGNATURES  JUDGED  EQUALLY  LOUD.  SOLID  LINE; 
LOUDNESS  OF  LOW-BOOM  SIGNATURES  VS  FLAT  TOP  DURATION. 


EXAMPLE  TIME-PRESSURE  HISTORY  OF  QUARRY  BLAST.  SIGNATURE 
WAS  RECORDED  AT  MILTON  LIMESTONE  QUARRY,  MAY  l6,  1974 
(from  REF.  U)  . SCALE;  1.78  psf /division,  125  ms/division. 


FIG.  29  REPRODUCTION  OF  "IDEALIZED"  BLAST  WAVEFORJG  BY  UTIAS  SIMULATION 
BOOTH  AS  USED  IN  THE  SUBJECTIVE  TEST.  RISE  TIME  = 0.22  ms, 
OVERPRESSURE  = 1 psf  (48  N/n^) . 

One  pulse,  duration  = 25  ms 
B)  Four  pulses,  duration  = 100  ms 
c)  Sixteen  pulses,  diiratlon  = 400  ms 
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FIG.  30  TRADE- OFF  BETWEEN  EFFECTIVE  LOUDNESS  AND  DURATION  OF  BLAST 
SIGNATURE.  COMPARISON  BETWEEN: 

1.  EXPERIMENTAL  RESULTS 

2.  CALCULATED  LOUDNESS  BY  JONNSON-RCBINSON  MEmOD  FOR 
Tc  = 175  ms 

3.  CALCULATED  LOUDNESS  BY  JOHNSON-RCBINSON  METHOD  FOR 
tq  = 70  ins 
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FIG.  31  COMPARISON  BETWEEN  EXPERIMENTAL  AND  CALCULATED  (JXXINSON-ROBINSON 

METHOD)  LOUDNESS  LEVELS  FOR  THE  HALF  N-WAVE  SIGNATURES  (CONTINUOUS 
LINE) . IDEAL  CORRELATION  BETWEEN  CALCULATED  AND  SUBJECTIVE  RESULTS 
IS  SHOWN  BY  BROKEN  LINE  (FROM  REF.  63) . 
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FIG.  32  TRADE-OFF  BETWEEN  EFFECTIVE  LOUDNESS  AND  DURATION  OF  BLAST  SIGNATURE. 
COMPARISON  BETWEEN; 

1.  E3CFERIMENTAL  RESULTS  (SOLID  LINE) 

2.  CALCULATED  LOUDNESS  BY  JOHNSON-ROBINSON  METHOD  WITH  BAND  PRESSURE 
LEVELS  DEFINED  AS 

(I/tjT)  r |F(a.)  I^dw  X (T/tc,  T < Tq;  1,  T > Tq)  (tq  = 175  ms) 

(BROKEN  LINE) 

3.  CALCULATED  LOUDNESS  BY  GENERALIZED  JOHNSON-ROBINSON  METHOD  WITH  BAND 
PELESSURE  LEVELS  DEFINED  AS 

U2 

(i/tiT)  r (F(u))I^du  X[l  . (t^  = 125  ms)  (BROKEN  LINE 
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FIG.  33  TRADE-OFF  BETWEEN  EFFECTIVE  LOUDNESS  AND  DURATION  OF  BLAST  SIGNATURE. 
COMPARISON  BETWEEN: 

1.  EXPERIMENTAL  RESULTS  (SOLID  LINE) 

2.  CALCULATED  LOUDNESS  BY  JOHNSON-RCBINSON  METHOD  WITH  BAND  PRESSURE 
mVEIfi  DEFINED  AS 


(I/tiT)  r X (t/tc  T < tc;  1,  T > Tc)  (tc  = 125  ms) 

(BROKEN  LINE) 


3-  CALCULATED  LOUDNESS  BY  GENERALIZED  JOHNSON-ROBINSON  METHOD  WITH  BAND 
PRESSURE  IJIVET^?  DEFTNEn  AS 


